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Abstract: This is a written account of the computer tutorial offered at the Sixth MC4BSM 
workshop at Cornell University, March 22-24, 2012. The tools covered during the tutorial 
include: FeynRules, LanHEP, MadGraph, CalcHEP, Pythia 8, HERWIG++, and 
Sherpa. In the tutorial, we specify a simple extension of the Standard Model, at the 
level of a Lagrangian. The software tools are then used to automatically generate a set of 
Feynman rules, compute the invariant matrix element for a sample process, and generate 
both parton-level and fully hadronized/showered Monte Carlo event samples. The tutorial 
is designed to be self-paced, and detailed instructions for all steps are included in this 
write-up. Installation instructions for each tool on a variety of popular platforms are also 
provided. 
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Figure 1. Simulation path from the Lagrangian £ of a BSM to a particle-level Monte Carlo event 
sample. All tools shown here are included in this tutorial. 

1 Introduction 

With the advent of the LHC, the field of particle physics has entered an exciting new era. A 
direct exploration of the TeV energy scale has finally begun. A large number of extensions 
of the Standard Model (SM), or "Beyond-the-Standard" Models (BSMs) for short, have 
been proposed over the years, and many of them can be probed by the LHC experiments. 
To make such tests possible, theoretical predictions of each model must be computed at a 
level that allows direct comparison to data. Monte Carlo (MC) generators are the basic 
tool for obtaining such predictions. 

Any BSM consistent with basic postulates of quantum mechanics and special relativity 
is a quantum field theory (QFT) defined by specifying its field content and the interaction 
Lagrangian C. As long as all interactions can be treated perturbatively, the fully differential 
cross section for any scattering process of interest can be computed by following the standard 
procedure of QFT perturbation theory. Namely, starting from C, one derives a set of 
Feynman rules, draws all possible Feynman diagrams for the process of interest, writes 
down the invariant matrix element Ai, squares it, and multiples it by the phase space 
weight. Since this procedure is completely unambiguous, it should be possible to teach 
a computer to perform it automatically. A set of automated tools to achieve this task 
has recently been developed, see Fig. 1. Feyn RULES [1] and LanHEP [2] perform the 
derivation of a complete set of Feynman rules from a given Lagrangian. The output of 
these programs can then be communicated to an automated matrix element calculator such 
as MadGraph [3] or CalcHEP [4] (or its sister package COMPHEP [5, 6]) which generates 
|A^P for any process specified by the user. This object in turns serves as an input to the 
Monte Carlo simulation of the process, first at the parton level with parton-level Monte 
Carlo event generators like MadGraph or CalcHEP (CompHEP). They create a set 
of "events", where each event is a record containing the identity and four-momentum of 
each of the initial- and final-state particles in the hard scattering process. The events are 
distributed in momenta, helicities, etc. according to the invariant matrix element 
computed in the previous step. In effect, the MC generator performs a pseudo-experiment: 
the generated set of events is just a particular statistical realization of the distributions 
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Table 1. The BSM field content (with quantum numbers) of the reference toy model. 

predicted by the theory for the process of interest. The produced parton-level events are 
in turn handed to a general purpose event generator such as Pythia [7] or Herwig [8, 9], 
which creates complete events, including the effects of fragmentation and hadronization 
of colored particles, initial and final state radiation via parton showers, effects from the 
underlying event, decays of unstable resonances, etc. The communication between the 
two classes of generators (parton-level and general purpose) is done following a universally 
accepted "Les Houches Accord" (LHA) format [10, 11]. 

The program of the Sixth MC4BSM workshop held at Cornell University in March 2012 
centered around hands-on computer tutorials [12] illustrating the tools depicted in Fig. 1. 
The starting point of the tutorials was a simple toy theory model described in Section 2. 
Its Feynman rules can then be automatically derived via either FeynRules (Section 3) or 
LanHEP (Section 4), the end product being the input files necessary to define the model in 
the parton-level event generators MadGraph and CalcHEP. The next step is to produce 
parton-level events, which is done in the tutorials covered in Sections 5 and 6, respectively. 
The parton-level events are then fed into Pythia 8 (Section 7) or HerwigH — h (Section 8). 
Finally, Section 9 contains a special tutorial on Sherpa [13], which provides an alternative 
path combining all of the above steps. In principle, each tutorial exercise is self-contained 
and independent of the others. Typically, each tutorial has two parts: in the first part 
(the pre-workshop exercise) students download and install the software and perform some 
simple tests to make sure it runs properly, while the second part deals with the actual 
physics simulations. 

2 The reference BSM model used in the tutorials 

The tutorial exercises are illustrated with a toy reference BSM model whose particle content 
is shown in Table 1. The model contains two real scalar fields, 0^ and 0^. They are singlets 
under all SM gauge groups. Their mass terms are^: 

2 2 

/^s.m. = -^0?-^02-"^?20102. (2.1) 

The corresponding mass eigenstates will be denoted by and <I>2, and their mass eigen- 
values by Ml and M2, respectively. For definiteness we will assume that Mi < M2. 

^AU Lagrangian parameters, here and below, are assumed to be real. 
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The model also contains two new Dirac fermion fields, U and E. Their SM quantum 
numbers are those of the SM ur and e/j, respectively. These fields have mass terms 

= MutJU + MeEE . (2.2) 

and interact with the new scalars via 

>CYuk = Xi^iUPru + \2<p2UPRU + A^i^P^e + >^2<p2EPRe , (2.3) 

where u and e are the SM up-quark and electron fields. Note that there is a Z2 symmetry 
under which all fields we added (</>i,2i U, E) flip sign, while all SM fields do not, so the new 
particles must be pair-produced, and the lightest new particle (LNP) is stable. This same 
Z2 also forbids U — u and E — e mixing via Yukawa couplings with the SM Higgs. 
We assume the following ordering of masses: 

Mu > M2 > Ml > Ml , (2.4) 

so that $1 is the LNP. Not having any SM interactions, it appears as MET in the detector. 
The ultimate goal of the tutorial is to simulate the process 

pp^UU, (2.5) 

at an 8 TeV LHC, and the subsequent U decays: 

[/ ^ u$i , (2.6) 
U^u^2, '^2^eE, E^e^i. (2.7) 



3 FeynRules tutorial 

The FeynRules tutorial consists of two steps: a pre- workshop exercise described in Section 
3.1, and the on-site exercise described in Section 3.2 

3.1 Installation instructions 

The aim of this session is to make the user familiar with the basics of FeynRules using 
the example of the Standard Model included in the distribution. FeynRules can be 
downloaded from http://feynrules.irinp.ucl.ac.be. You need to have Mathematica 
installed on your machine in order to run FeynRules. Simply download the package and 
untar it. 

Next, you can have a look at the SM implementation, contained in the directory 
/f eynrules-current/models/SM/. This folder contains various files, but in this intro- 
ductory session we will only look at two of them: 

1. SM.fr: the model file (a text file) containing the implementation of the SM in Feyn- 
Rules. 

2. SM.nb: a Mathematica notebook, showing how to load and run the model to obtain 
the Feynman rules. 
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3.1.1 The FeynRules model file 

Open the model file for the SM, SM.fr, in a text editor. In the following we give a very 
brief account on the general structure of a FeynRules model file. During the tutorial, you 
will be given the opportunity to write your own FeynRules model file, and we therefore 
limit ourselves at this stage to describe the generic features. 
A generic FeynRules model file contains four parts: 

1 . A list called M$GaugeGroups containing the definitions of all the gauge groups of the 
model. 

2. A list called M$Parameters containing the definitions of all the parameters of the 
model. 

3. A list called M$ClassesDeclarations containing the definitions of all the fields of the 
model. Note that the fields are grouped into classes of particles carrying the same 
quantum numbers. 

4. The Lagrangian of the model, written in Mathematica form. 

3.1.2 The sample notebook 

Next, open the notebook SM.nb. The first few lines of the notebook show how to load 
FeynRules into Mathematica, 

$FeynRulesPath = SetDirectory ["< address of the package >"] 
« FeynRules' 

SetDirectory [$FeynRulesPath <> "/Models/SM"] ; 

The first line sets the path to the FeynRules main directory. Note that < address of the package > 
should be replaced by the path of the FeynRules main directory on your computer. For 
example, if the FeynRules main directory is in your home folder, then the first line should 
read 

$FeynRulesPath = SetDirectory ["~/feynrules-current"] 

The second line loads the FeynRules package. If successful, an output is printed on the 
screen. The third line changes the path to the subfolder Models/SM that contains the SM 
implementation. 

After FeynRules has been successfully loaded into Mathematica, the SM implemen- 
tation can be loaded by issuing the command 

LoadModel ["SM.fr"] 

The next command is optional, and loads some additional files that put the masses of the 
first two generations of fermions to zero, and restricts the CKM matrix to be diagonal. As 
this step is optional, we will not describe it any further. 
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The first section of tlie notebook, The SM Lagrangian, sliows tlie Lagrangian tliat is 
implemented in SM.fr. Note tliat tlie SM implementation of FeynRules is available in 
both Feynman and unitary gauge. The boolean variable FeynmanGauge allows to switch 
between the two gauges. The second section shows how one can use some built-in functions 
to perform some basic checks on the Lagrangian, e.g., to check if the Lagrangian is hermi- 
tian. The last two sections illustrate how to compute the Feynman rules of the model with 
FeynRules, and how to use the interfaces to the Feynman diagram generators. 

3.2 On-site tutorial 

The aim of this tutorial is to implement a simple extension of the SM into FeynRules. 
The model itself is described in Section 2. In the following we give detailed instructions of 
how to implement this model into FeynRules and how to obtain the corresponding model 
files for CalcHEP, MadGraph 5 and Sherpa. 

3.2.1 Preparation of the model file 

As the model we are going to implement is a simple extension of the SM, it is not necessary 
to start from scratch, but we can use the implementation of the SM included in the folder 
/Models/SM. We therefore start by making a copy of this folder, in order to keep a clean 
version of the SM. To do so, change directory to the Models subdirectory and make a copy 
of the SM folder, before going into the new directory 

cd Models 

cp -r SM MC4BSM 

cd MC4BSM 

Even though the implementation is based on the model file SM.fr for the Standard Model, 
the SM sector of the model will be implemented into a separate file that will simply be 
loaded on top of SM.fr. We therefore start by opening a blank text file called MC4BSM.fr. 
You can start by personalizing the model file by including a name for you model, the name 
of the author, etc., 

M$ModelName = "MC4BSM_2012" ; 

M$Iiif ormation = {Authors 
Version 
Date 

Institutions 
Emails 

}; 

Note that his information is purely optional and could be omitted. 

3.2.2 Implementation of the new parameters 

We start by implementing the new parameters. The model we are considering depends on 
9 new parameters, which we assume to be real in the following. FeynRules distinguishes 



-> {"C. Duhr"}, 

-> "1.0", 

-> "27. 02. 2012", 

-> {"ETH Zurich"}, 

-> {"duhrc@itp.phys.ethz.ch"} 
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between two types of parameters, the so-called external parameters given as numerical 
inputs and the internal parameters, related to other external and/or internal parameters 
via algebraic expressions. All the parameters of the model, both external and internal, are 
given in the Feyn Rules model file as a list named M$Parameters. Note that if an internal 
parameter x depends on some other parameter y, then y should appear in M$Parameters 
before the internal parameter x. 

The new external parameters of the model are 

• 5 mass parameters: mi, m2, mi2, Mjj, Me- 

• 4 coupling constants: Ai, A2, A'^, A2. 

Note however that there is a difference between the mass parameters in the scalar and 
fermonic sectors: while the masses in fermonic sector are physical masses, the mass matrix 
in the scalar sector is not diagonal. For this reason, we will not discuss in the following 
the fermion masses Mjj and Me, as they will be defined together with the particles rather 
than as parameters of the model. 

Let us now turn to the definition of the mass parameters in the scalar sector. The 
masses mi, m2 and mi2 will be denoted in the Feyn RULES model file by MMl, MM2 and 
MM12. In the following we only show how to implement MMl, all other cases being similar. 
MMl corresponds to the following entry in the list M$Parameters, 

M$Parameters = { 

MMl == { 

ParameterType -> External, 
Value -> 200 

}. 

} 

The first option tags MMl as an external parameter, while the second option assign a value 
of 200GeV to mi . We stress that this numerical value can be changed later on in the matrix 
element generators. 

The masses in the scalar sector are not the physical masses, because the mass matrix 
is not diagonal. In order to obtain the physical masses, we need to diagonalize the mass 
matrix 

fi^iy (3.1) 

In the following, we denote the eigenvalues by MPel and MPe2. In addition, we need to 
introduce a mixing angle 9 (th) relating the fields (f)i to the mass eigenstates $i by. 
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As in this case the mass matrix is only two-dimensional, we can compute the eigenval- 
ues and the mixing angle analytically, and simply implement the analytical formulas into 
FeynRules. The implementation follows exactly the same lines as for the masses mi, m2, 
mi2, with the only differences that 

1. the ParameterType is Internal (as these parameters are dependent on the external 
mass parameters, 

2. the Value is given by an analytical expression (in Mathematica syntax). 

Next we turn to the implementation of the new coupling constants, which we will call 
laml, la21, lamlp, lam2p. They are all external parameters, and thus the implementation 
follows exactly the same lines as the implementation of the mass parameters, with only one 
modification: some matrix element generators, like for example MadGraph, keep track of 
the types of couplings that enter a process. This allows for example to generate a process 
by only taking into account QCD-type vertices, and to neglect all QED-type vertices. For 
this reason, it is mandatory to tell the matrix element generator how the new coupling 
constants should be counted. As in this case we are dealing with new classes of couplings 
which are a priori independent of QCD or QED interactions, we simply assign a new tag, 
called interaction order, to the coupling via the option 

Interact ionOrder -> {NP, 1} 

The name of the tag (NP for "new physics" in this case) can be chosen freely. The above 
option instructs the matrix element generator to count one unit of "NP" for each new cou- 
pling. 

3.2.3 Implementation of the fields 

In this section we discuss the implementation of the new fields. The implementation is 
similar to the implementation of the parameters, i.e., all the fields are entries of a list 
called M$ClassesDescription. In Tab. 2 we show the names of the fields used in the 
implement at ion ^ . 
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01 


(t>2 
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$2 


uv 


ev 


pil 


pi2 


Pl 


p2 



Table 2. Symbols used for the fields in the FeynRules implementation. 

We illustrate the implementation of a new field on the example of the particle U (uv). 
The definition of the particle corresponds to an entry in M$ClassesDescription of the 
following form 

^Note that the symbol u, e and phi are aheady in use in the SM implementation. We also avoid using 
simply uppercase letters, as some matrix element generators are case insensitive. 
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M$ClassesDescription = { 



F[10] == { 



ClassName 



> 



uv, 



Self Conjugate 
Indices 



> 



> 



False , 

{Index [Colour] } , 
{Y -> 2/3, Q -> 2/3}, 
{Muv, 500}, 
{Wuv,l} 



QuantumN umbers 



> 



Mass 



> 



Width 



> 



}. 



} 



The meaning of this definition is as follows: each particle class has a name of the form 
X[i], where X is related to the spin of the field (See Tab. 3), and i is an integer that 
labels the classes. Note that i can be chosen freely, as long as there is no name clash with 
an already existing class (in this case, there could be a name clash with the SM particles 
already defined in SM.fr). Each class has a series of options 

1. ClassName: the symbol by which the particle will be represented in the Lagrangian. 

2. Self Conjugate: a boolean variable, indicating whether the particle has an antiparti- 
cle (False) or not (True). If the field is not self-conjugate, a symbol for the antiparticle 
is automatically defined by appending "bar" to the name of the particle. In the above 
example the antiparticle associated to uv will be denoted by uvbar. Note that in the 
case of fermions the symbol for the antiparticle refers to the quantity U rather than 



3. Indices: All indices carried by the field. The available types of indices from the SM 
implementation are 

• Generation: fermion flavor index ranging from 1 to 3, 

• Colour: fundamental color index ranging from 1 to 3, 

• Gluon: adjoint color index ranging from 1 to 8, 

• SU2W: adjoint SU{2)l index ranging from 1 to 3. 

4. QuantumNumbers: a list of all U{\) charges carried by the field. In the SM implemen- 
tation the following U{1) charges are already defined 

• Y: weak hypercharge, 

• Q: electric charge. 

5. Mass: the mass of the particle. It is a list of two elements, the first being the symbol 
used to represent the mass, and the second its value (in GeV). If the value of the mass 



C/t. 
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is obtained from some analytic expression defined as an internal parameter with the 
same symbol (as is the case for example in the scalar sector of the model), the value 
is set to Internal. 

6. Width: the width of the particle. The definition is similar to Mass. Note that as we 
do not yet know the widths of the new particles, we simply set it for now to IGeV, 
and will determine its exact value later using one of the matrix element generators. 

The implementation of the other mass eigenstates (ev, pi, p2) is similar, so we do not 
discuss it here. 



Spin 


1/2 1 2 ghost 


Symbol 


S F V T U 



Table 3. Available particle classes in FeynRules. 



Let us comment on the implementation of the interaction eigenstates Indeed, while 
the matrix element generators work exclusively at the level of the mass eigenstates, the 
interaction eigenstates are in general useful to write the Lagrangian in a compact form. It 
is therefore useful to define also the fields for the interaction eigenstates (pi. The definition 
of these fields is similar to the mass eigenstates, e.g., 

S[10] == { 

ClassName -> pil, 
Self Conjugate -> True, 
Indices -> {}, 

Unphysical -> True, 

Definitions -> {pil -> - Sin[th] pi + Cos [th] p2} 
}. 

First, note that the Mass and Width options are omitted'^, as these fields are not mass 
eigenstates. This last fact is made explicit by the option 

Unphysical -> True, 

which instruct FeynRules not to output this field to a matrix element generator. Finally, 
the relation of the field pil to the mass eigenstates is simply given as a Mathematica 
replacement rule in the Definitions option. 

3.2.4 Implementation of the Lagrangian 

The definitions in the model file being complete, we now turn to the implementation of 
the Lagrangian. This can be done either in the model file, or directly in a Mathematica 
notebook. Here we use the latter approach, and we start by opening a new notebook and 
load the FeynRules package (see the pre-installation instructions). Next we have to load 
the model files, both for the SM and for the new sector, 



^The QuantumNumbers option is also omitted, but for the simple reason that the fields cpi do not carry 
any (7(1) charges. 
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LoadModel [ " SM . f r " , "MC4BSM . f r "] 

Note that the new model file should be loaded after SM.fr. Furthermore, we also load 
two additional files, which restrict the first two fermion generations to be massless and the 
CKM matrix to be diagonal, 

LoadRestriction["DiagonalCKM.rst" , "Massless. rst"] 

The new Lagrangian consists of three parts, 

~ ^scalar,kin ~l~ ^ fermion,kin "I" ^Yuk • ("^•"^) 

The kinetic terms for the new scalars can easily be implemented by using the symbols for 
the gauge eigenstates and the mass parameters defined in the model file, as well as the 
symbol for the space-time derivative dfj, in FeynRules, del [ . . . , mu] . As an example, 
we have 

1/2 del [pi, mu] del [pi, mu] - 1/2 MMl-2 pl-2 

The kinetic terms for the fermions can be implemented in a similar way. However, as the 
fermions are charged under the SM gauge group, we have to use the covariant derivative 
DC rather than the space-time derivative del. Furthermore, we have to use a "." instead of 
an ordinary multiplication in order to take the non-commuting nature of the fermions into 
account. As an example, we have 

i U-i^D^U - MutJU 
I uvbar . Ga [mu] . DC [uv , mu] - Muv uvbar . uv 

where Ga[mu] is the FeynRules symbol for the Dirac matrix 7^^. Finally, the Yukawa 
interactions can be implemented in the same way as the kinetic terms for the fermions, e.g., 

Ai tJP+u 
laml pil uvbar . Pro j P . u 

where u denotes the u quark field defined in SM.fr and ProjP denotes the right chiral pro- 
jector (the left projector is denoted by ProjM). Note that FeynRules contains a function 
HC [ ] which allows to obtain the hermitian conjugate of an expression in an automated 
way. 

3.2.5 Computing the Feynman rules and running the interfaces 

Our model implementation is now complete, and so we can compute the Feynman rules. 
The Feynman rules of the new sector can be obtained by issuing the command 

FeynmanRules [ LNew ] 

where LNew is the name of the variable that contains the new Lagrangian. 

The Feynman rules can be written to file in a format suitable to various matrix element 
generators by using the FeynRules interfaces. In this tutorial, we will use the interfaces 
to CalcHEP, MadGraph 5 and Sherpa, which can be called via 
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WriteCHOutput [ LSM + LNew ] ; 
WriteUFO [ LSM + LNew ] ; 
WriteSHOutput [ LSM + LNew ] ; 

where LSM is the SM Lagrangian implemented in SM.fr. Note that the SM implementation 
is available in both Feynman gauge and unitary gauge. A boolean variable FeynmanGauge 
allows to switch between both gauges. While CalcHEP and MadGraph 5 support both 
gauges, Sherpa only supports unitary gauge, and the FeynmanGauge variable should cor- 
respondingly be set to False. 

4 LanHEP tutorial 

LanHEP is a program which allows to create input model files for the programs COMPHEP 
[14] or CalcHEP [15] from a Lagrangian provided by the user. This program is as easy as 
M^]X. We define the symmetries, the particles and then write down a Lagrangian, just as 
we would write them in a paper. Then we compile the write-up and obtain COMPHEP- or 
CALCHEP-specific model files for further collider studies. In short, if you can write down 
a Lagrangian, and if you have enough knowledge of IM^i]X, then your skills are probably 
enough to use LanHEP as well. 

An excellent reference for beginners is the LanHEP manual [16]. Here we will demon- 
strate how to use LanHEP with the specific example (see Section 2) provided for the 
MC4BSM 2012 workshop''. First let us discuss how to install LanHEP. 

4.1 Installation 

Installation is never an easy job, especially when there are several choices for the operating 
system. But if you have a working C-compiler in your operating system, there should be 
no problem. 

1 Linux with gcc. 

2 OS/X: to use gcc, you need to install Xcode in Apple Appstore. But to use Monte 
Carlo simulations, you will also need a fortran compiler. Free fortran compilers are 
available from these sites [18, 19].^ 

3 Windows: we recommend that you use a virtual machine^ to install Linux under the 
Windows system or have a dual boot with a Linux partition. Most Monte Carlo 
simulations are too complicated to run under Cygwin. 

The detailed installation instructions are well explained on the LanHEP webpage [17]. 
Here we will briefly recap the procedures based on the file, "lanhep315.tar" which is the 
current version at the time of this writing. 

*MC4BSM 2012 provides a nice homepage [12] for all tutorials. 

^Starting from Mountain Lion, in Xcode, you need to install "Command Line Tools" which you can find 
in "Xcode/Preferences/Downloads". The current version coming from Xcode is gcc 4.2, thus we recommend 
that you match your versions of gcc, gH — h and gcov with the version of gfortran that you download from a 
third party. 

®There is a free virtual machine software called VirtualBox [20]. 
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1. Download the LanHEP program from the LanHEP webpage [17]. 



2. Make a proper dh'ectory for the LanHEP installation. For example 

mkdir $HQME/MC4BSM 

3. Move the downloaded file to the directory that you just created: 

a. Let's suppose that you made an MC4BSM directory under your $HOME di- 
rectory, and you downloaded a "lhep315.tar" file in the '$HOME/Downloads" 
directory, then you can move the file into your MC4BSM folder by typing 

mv $H0ME/Downloads/lhep315.tar $HQME/MC4BSM/ . 

b. Unzip "lhep315.tar" file by typing 
tar -xvf lhep315.tar 

3. Go to the "lanhep315" directory and read the "README" file. 

4. Compile LanHEP by typing 

make 

4.2 Warm-Up exercise 

Before we begin the real MC4BSM exercise, it would be good to start with a very simple 
example, "QED", which is the first example in the LanHEP manual version 3.0 [2]. You 
can find its model files in the "mdl" directory as "qed.mdl" . The contents of this file are 
the following 

model qed/ 1 . 

parameter ee = 0.3133: 'Electric charge', 
vector A/A: photon. 

let F~mu"nu=deriv~mu*A"nu-der iv~nu*A~mu . 
spinor el : (electron, mass me=0 . 000511) . 
Iterm ee*El*gamma*A*el . 

Without looking at the manual, we can easily guess that there are two fields, an electron 
(el) and a photon (A). In this particular model file, there are missing terms from the 
Lagrangian (if you compare it to the text in the manual). You can modify this model file 
by writing down the missing terms. Thus a complete simple "QED" model file will look like 

model qed/ 1 . 

parameter ee=0. 31333: 'elementary electric charge', 
spinor el/El : (electron, mass me=0 . 000511) . 
vector A/A: (photon) . 

let F~mu"nu=deriv~mu*A"nu-der iv~nu*A~mu . 

Iterm -1/4* (F~mu~nu) **2 - 1/2* (deriv~mu*A~mu) **2 . 

Iterm El* (i*gairana*deriv+me) *el . 

Iterm ee*El*gairana*A*el . 

At this stage we can ask some simple questions: 
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1. What is the relation between el and El? 

2. How can we use a covariant derivative instead of sphtting the Lagrangian into pieces?''' 

3. What is the unit of mass in LanHEP? 

In order to create the model files for the actual COMPHEP or CalcHEP program, we 
need to compile this model file. It would be better to make a separate directory for the 
"compiled" model. Let's make the directory inside your LanHEP directory (if you are in 
the mdl directory, go one level up). For example, 

mkdir QED 

Now compile the "qed.mdl" file by using the "Ihep" command as follows^ 
./Ihep -QutDir QED mdl/qed.mdl -CalcHEP 

You will get this "warning" message 
File mdl/qed.mdl processed, sec. 

Warning: property 'pdg' is not defined for particle 'el' 

"PDG" is a number which is defined for each particle according to the Particle Data Group 
for various Monte Carlo simulators. For example the PDG number for an electron is 11. 
Our particle "el" defined in the "qed.mdl" file was not assigned a PDG number. There are 
two ways to resolve this warning. If you open the "calchep.rc" file in the mdl directory, you 
will see the following contents. 

prtcformat f ullname : 'Full Name ', 

name : ' P ' , aname : ' aP ' , pdg : ' number ' , 

spin2, mass, width, color, aux, texname: ' LaTeX(A) ', 

atexname:' LateX(A+) '. 

prtcproperty pdg: (A=22, Z=23, 'W+'=24, G=21, 

d=l, u=2, s=3, c=4, b=5, t=6, 
ne=12, nm=14, nl=16, 
e=ll, m=13, 1=15, 
~ne=1000012, ~nm=1000014, ~nl=1000016, 
~el=1000011, ~ml=1000013, ~11=1000015, 
~e2=2000011, ~m2=2000013, ~12=2000015, 
~eL=1000011, ~mL=1000013, 
~eR=2000011, ~mR=2000013, 



~ul=1000002, ~cl=1000004, ~tl=1000006, 

~u2=2000002, ~c2=2000004, ~t2=2000006, 

~uL=1000002, ~cL=1000004, 

~uR=2000002, ~cR=2000004, 



^See Section 3.4 in the manual 3.0 [2]. 

*If you type "./Ihep" without any options, you will see the available choices. 
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~dl=1000001, ~sl=1000003, ~bl=1000005, 

~d2=2000001, ~s2=2000003, ~b2=2000005, 

~dL=1000001, ~sL=1000003, 

~dR=2000001, ~sR=2000003, 

h=25, H=35, H3=36, 'H+'=37, 

~ol=1000022, ~o2=1000023, ~o3=1000025, ~o4=1000035, 
'~1+'=1000024, '~2+'=1000037, ~g=1000021) . 

Thus, we can edit this file directly to include a PDG number (11) for our field 'el', or we 
can use a function "prtcproperty pdg" during our write-up of the Lagrangian. 

4.2.1 Output 

After compiling the model file "qed.mdl", we will get four files in the QED directory. 

1. Igrngl.mdl : This file contains the Feynman rules. 

2. prtclsl.mdl: This file contains the list of particles and their properties, for example, 
spin, mass and width. 

3. varsl.mdl: This file contains parameters that we defined, for example, the gauge 
coupling "ee" and the mass of the electron "me". 

4. funcl.mdl: If we want to introduce a function which depends on the above parameters, 
then this file contains the definition of this function. 

Since these files are specific to CompHEP and CalcHEP, it would be good to consult 
their manual [4] and the TASI 2012 lecture note [21] for more detailed information. 

4.3 MC4BSM 2012 on-site exercise 

In the MC4BSM 2012 workshop tutorial, we work with the simple model described in 
Section 2. In this section, we will show the additional Lagrangian terms for this model. 



with corresponding kinetic terms. LanHEP already contains a model file for the SM in the 
"mdl" directory, thus it is convenient to start from this pre-existing model file, in our case, 
"newsm.mdl" file. To write down a model file using LanHEP, we need to understand 

1. How to read a pre-defined model file and how to use this in our model file. 

2. How to relate mass eigenstates to interaction eigenstates. 

3. How to check whether our model file is written up correctly or not. 




■s.m 



'"l ,2 '"'2 i2 2 I I 
- -^'^2 - "ll2'Pl<?'2, 

MuUU + MeEE, 

XkPiUPru + X2(t)2UPRU + X^c^iEPrb + X'2hEPRe + h.c. 



(4.1) 
(4.2) 
(4.3) 
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4.3.1 Setup 



Let's name our file "mc4bsm.mdl" and save it in the "mdl" directory since this directory 
contains files that we will be using. 

read newsm. 

keys CKMdim=l. 

keys SMmassless=l . 

keys gauge_f ixing=Feyiiman. 

y, Remove the 1st and 2nd generation mass 
do_if SMmassless==l . 
parameter Mm = , Ms =0, Mc=0 . 
end_if . 

0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/ 

/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/^ 

yoyoy.yoyo7.yoyo7oyo7.yoyo7oyo7.y.yo7oyo7.y.yoyoyo7.yoyoy.yoyoyoyoy.yoyoyoyo7.7o 

7. Make CKM a diagonal Matrix. 
do_if CKMdim==l. 

parameter Vub=0, Vcb=0, Vtd=0, Vts=0, Vtb=l, Vud=l, Vus=0, Vcs=l, Vcd=0. 
end_if . 

0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/^ 

/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/^ 

do_if gauge_f ixing==Feynman. 

model 'MC4BSM_f eynman'/4. 
do_else_if gauge_f ixing==unitary . 

model 'MC4BSM_unitary'/3. 
end_if . 



On the first line in this code, we see the "read" command. This command allows us to 
use a model file which is already defined in the "mdl" directory, in our case "newsm. mdl". 
There are also lines starting with "keys" and below them a section with "do_if, end_if". 
By using this method, we can control different set-ups for our model file. If you look into 
"newsm. mdl" file, the CKM entries. Gauge choice and masses of SM particles are already 
defined. But after we read this model file, we can over- write/replace those pre-defined 
parameters with new values of our choice. We put the CKMdim and SMmassless flags since 
the FeynRules example from Section 3 uses restrictions on the CKM and Particle masses. 
In above code, we can also see that "%" is a flag to comment out a line. 

4.3.2 Particles and corresponding parameters 

We will follow the same conventions as the FeynRules tutorial in Section 3 (see Table 4). 



Fields 


U 


E 




4>2 


$1 


<I>2 


LanHEP 


uv 


ev 


pil 


pi2 


Pl 


p2 



Table 4. Symbols used for the fields in the LanHEP implementation. 
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For the scalar fields, (<I>i, <I>2) are the mass eigenstates (physical states) of {(pi, (j)2) with 
mixing angle 6: 



— sin cos 9 
cos 9 sin 9 



$1 

$2 



(4.4) 



coming from the eigen states of the following mass matrix, 



mT m 



12 



m\2 m-2 



(4.5) 



Thus there are five mass parameters, including the two eigenvalues (MPel, MPe2) of 
the mass matrix, and a corresponding mixing angle 9, summarized in Table 5. 



Mass parameters 


mi 


1712 


mi2 


Mu 


Me 


tan((9) 


LanHEP 


MMl 


MM2 


MM12 


Muv 


Mev 


tth 



Table 5. Symbols used for the mass parameters in the LanHEP implementation. 



For coupling constants, we use the notation in Table 6. 



Coupling constants 


Ai 


A2 


a; 


A'2 


LanHEP 


laml 


lam2 


lamlp 


lam2p 



Table 6. Symbols used for coupling parameters in the LanHEP implementation. 



We sum up the above set-up for the particles in the following: 
'/, Define new parameters. 

parameter MMl = 200 : 'Mass parameter for pil', 
MM2 = 300 : 'Mass parameter for pi2', 
MM12 = 50 : 'Mixing term for pil and pi2', 
laml = 1 : 'Yukawa interaction for pil, U and u' , 
lam2 = 1 : 'Yukawa interaction for pi2, U and u' , 
lamlp = 1 : 'Yukawa interaction for pil, E and e', 
lam2p = 1 : 'Yukawa interaction for pi2, E and e', 
Muv = 500.0 : 'Mass parameter for uv and v' , 
Mev = 250.0 : 'Mass parameter for ev and e'. 

parameter 

MPel = sqrt (MM2**2+MMl**2-sqrt (4 . 0*MM12**4+ (MM2**2-MM1**2) **2) ) /Sqrt2 , 
MPe2 = sqrt (MM2**2+MMl**2+sqrt (4 . 0*MM12**4+ (MM2**2-MM1**2) **2) ) /Sqrt2 . 

parameter 

tth = (MM2**2-MMl**2+sqrt(4.0*MM12**4+(MM2**2-MMl**2)**2))/(2.0*MM12**2) . 
parameter sth = tth/sqrt (l+tth**2) , cth = 1 . 0/sqrt (l+tth**2) . 
angle sin=stli,cos=cth,tan=tth, texname='\\theta' . 

Define new particles, 
spinor ev: ('heavy electron' , mass Mev, width Wev=auto) , 
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uv: ('heavy quark', color c3, mass Muv, width Wuv=auto) . 

scalar pl/pl : ( 'LNP ' ,mass MPel) , 

p2/p2 :( 'Heavy scalar', mass MPe2, width Wpe2=auto) . 

let pil= -sth*pl+cth*p2, pi2=cth*pl+sth*p2 . 

prtcproperty pdg: ( ev=9000009,uv=9000008,pl=9000006,p2=9000007) . 

As we notice here, the mathematical commands are similar with Fortran, for example 
— )• X * *2. For more information, we refer the user to the manual v. 3.0 [2]. There are 
three major points that need to be explained: 

1. Particle width: CompHEP/CalcHEP can calculate the width of a particle on the fly. 

Thus in the particle definition, we can set it to "auto" and then CompHEP/CalcHEP 
will calculate it. 

width Wev=auto 

2. Physical state: We defined the physical states (^1,^2) in the particle definitions and 

can use the "let" command to define interaction eigenstates: 

let pil= -sth*pl+cth*p2, pi2=cth*pl+sth*p2. 

We will define our Lagrangian in terms of interaction eigenstates and LanHEP will 
automatically replace terms with the mass eigenstates. A Width of pi is not defined, 
because pi is a stable particle in our model. By default, its width = 0. 

3. PDG number: We set the PDG number to be consistent with the model from Feyn- 

RULES. In this onsite exercise, we will not change the "calchep.rc" file, but use the 
following command instead 

prtcproperty pdg:( ev=9000009,uv=9000008,pl==9000006,p2=9000007). 
4.3.3 Setting up the Lagrangian 

We start this subsection with the code for the Lagrangian, which we will then explain step 
by step. 

y. Kinematics terms for new fermion particles 

y, New fermion field is a counter part of SH right-handed particle 
let PR = (l+g5)/2. 

Iterm anti (psi) *gamma*PR* (i*deriv - Y*gl*Bl)*psi 
where 

psi=ev, Y= -1; 
psi=uv, Y= 2/3. 

Iterm GG*anti (uv) *lambda*gaimna*G*uv . 
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% Kinematic terms for scalars 

Iterm deriv*pil*deriv*pil/2+<ieriv*pi2*<ieriv*pi2/2. 
'/» New mass term 

Iterm -MMl**2/2 * pil*pil- MM2**2/2 * pi2*pi2- MM12**2*pil*pi2 . 
Iterm -Muv*anti (uv) *uv-Mev*anti (ev) *ev . 

% New Yiikawa type interaction. 

Iterm laml*pil*anti (uv) *PR*u+lam2*pi2*aiiti (uv) *PR*u 

+lamlp*pil*anti(ev)*PR*e+lam2p*pi2*anti(ev)*PR*e + AddHermConj . 

% Checking a La^angian 
CheckHerm. 



Fermions. 

1. We defined a right-handed projector PR with g5 (g5 is defined in newsm.mdl as g5=75.) 

2. We used 'Svhere" grammar to reduce unnecessary repeats and to give a clear structure. 

• psi is a generic notation for fermions here. Bl is U(1)y gauge field defined in 
"newsm.mdl" file. 

• anti(i/;) denotes ip, similarly to the FeynRules conventions. 

• Y is the hypercharge of the new particles, and gl is defined to be the U(1)y 
coupling in "newsm.mdl". 

• For uv(colored particle) there is an interaction with G (Gluon) through the 
strong coupling GO defined in "newsm.mdl" . 

For example, 

Lagrangian (9^ - 51 ^i//) V' 

LanHEP anti(psi)*gainma*PR*(i*deriv - Y*gl*Bl)*psi 



Notice that we can omit the Lorentz indexes when they are summed up. Alternatively, 
as in the "qed.mdl" file, we can write them down explicitly: 

anti(psi)*gamma"mu*PR*(i*deriv"mu - Y*gl*Bl~mu) *psi 

Scalar pgirticles. Usually, it is easy to write down the Lagrangian in terms of interac- 
tion eigenstates, instead of mass eigenstates (or physical states). Thus when we define 
Lagrangians, we use interaction eigenstates. During compiling of LanHEP, the interaction 
states will be replaced with physical states and we will not see unphysical states in the 
CompHEP/CalcHEP model files. 

To summarize, the procedure is the following: 

1. Define physical states. 
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2. Relate unphysical states to the corresponding physical states using the "let" command. 

3. Describe Lagrangian in terms of unphysical states. 

The same technique is applied in "newsm.mdl" for gauge fields, thus we recommend 
the user to read the "newsm.mdl" file at this stage. 

4.3.4 Finalizing the Lagrangian script 

We add 

+ AddHermConj 

to the Yukawa interaction to get a hermitian conjugate term. To check our script for the 
Lagrangian, we add the 

CheckHerm. 

command to make sure that our Lagrangian is hermitian. For additional commands for 
cross-checks, see section 6 of the LanHEP manual 2.0 [2]. 

We can get a model file for CompHEP/CalcHEP by compiling in the LanHEP 
directory 

mkdir tutorial 

./Ihep -OutDir tutorial mdl/mc4bsm.mdl -CalcHEP 

SO that we do not mix up our newly created model files with existing "mdl" files. 

5 MadGraph tutorial 
5.1 Preamble 

MadGraphS [22] can be run on a local computer or via the web at one of the following 
website: 

• http://madgraph.hep.uiuc.edu 

• http:/ /madgraph. phys.ucl.ac.be 

The registration is straightforward, and you can instantly creates optimized code for the 
computation of the cross-section of any processes. That code can be run directly on the 
web or downloaded and run locally. However, for security reason, generating events on the 
web is allowed only after that you have sent an email to one of the authors of MadGraph5. 
Since most of the functions are available on the internet, most users will not need to install 
MadGraphS. 

The MadGraphS collaboration plans to continue to improve this code both by adding 
new functionalities and by making the code easier to use. In particular the way to compute 
decay widths presented in this paper will soon be improved. New tutorials and updates on 
this tutorial can be found at [23]. 
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5.2 Installation of MadGraphS 



5.2.1 Installation on Windows 

MadGraphS and the associated programs are designed and tested on Linux and MacOs 
operating systems. The windows compatibility via cygwin is currently not supported. For 

Windows user, we advise they install Linux in dual boot. Virtual machines are another 
possibility. Note that some virtualbox packages do not include the library readline. This 
library is not mandatory but enables the auto-completion (See the paragraph associated to 
python installation to learn how to solve this). 

5.2.2 Installation on Linux / MacOS 

MadGraph 5 The last version of MadGraph5 can be found at the following page: 
https:/ /launchpad.net/madgraphS. This website is also the place, where you can ask ques- 
tion, make suggestions or report a bug. The installation is straightforward since you have 
only to untar it 

tar -xzpvf MadGraph5_vl .X.Y.tgz 

No compilation are required for MadGraphS, you can just launch it. 
. /MadGraph5_vl . X . Y/bin/mg5 

If you don't have a valid python version, MadGraph S will crash directly with an explicit 
message. In this case, you will need to install Python2.7. 

If you have admin rights on your system, you can run the following command: 

sudo In -s MadGraph5_vl .X. Y/bin/nig5 /usr/local/bin 
such that MadGraphS can be launched from any directory. 

Python The only requirement for MadGraphS is to have a current version of PYTHON (ver- 
sion 2.6 or 2.7). In most cases, it can be installed via your favorite repository manager. 
However, some of the linux repository distributes python 2.5. In that case, you will need 
to download python at the following link: http://www.python.org/download/ and follows 
the associate instructions. 

Note that for some linux versions (especially in virtual machine), the library readline 
is not present on the system. This will disable the auto-completion. If you care about that 
point, you will should first install that library (via your repository) and then to recompile 
python from the source code. 



- 21 - 



Optional package Various optional packages can be linked to MadGraph5 in order to 
customize the output, create plots, ... The installation of those packages is easy since you 
can install them by launching mg5 and typing 

nig5> install NAME 

Where NAME is one of the following package name: 

• MadAnalysis: A package to draw automatically various histogram linked to the event 
generation. 

• ExRootAnalysis [24]: A package to convert the various output in a ROOT format. 

• pythia-pgs: A package containing Pythia6 [25] and PGS [26]. Pythia6 is able to 
shower and to hadronize your events and is able to perform the matching [25, 27] for 
multi-jet production. PGS is a fast detector simulation package. 

• Delphes [28]: A package allowing to have a fast detector simulation, in replacement 
of PGS. 

Note that some of those programs might have some extra-dependencies (especially in Root) . 

Additional instructions for MacOs Compared to Linux the installation on MacOS 
might be more complex since MacOS doesn't provide various set of default programs present 
on Linux. Two important programs which are not present by default are gmake and gfor- 
tran4.x. We advise you to first check if those program are install via the commands. 

$> make --version 

$> gfortran --version 

In order to install gmake it is easiest to install xcode (free but requiring an apple developer 
account) 

• MacOs 10.5: https://connect.apple.com/cgi-bin/WebObjects/ 
MemberSite.woa/wa/getSoftware?bundleID=20414 

• MacOs 10.6: http://connect.apple.com/cgi-bin/WebObjects/ 
MemberSite.woa/wa/getSoftware?bundleID=20792 

• MacOs 10.7: http://itunes. apple.com/us/app /xcode/id448457090?mt=12 
Concerning gfortran you can download it from the following link: 

• MacOs 10.5: http://prdownloads.sourceforge.net/hpc/ 
gcc-lion . t ar . gz ?download 

• MacOs 10.6: http://prdownloads.sourceforge.net/hpc/ 
gcc-snwleo-int el-bin . t ar . gz ?download 

• MacOs 10.7: http://sourceforge.net/projects/hpc/files/hpc/gcc/ 
gcc-leopard-intel-bin.tar .gz / download 
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5.2.3 Testing the installation and learning the syntax 

MadGraphS includes a build-in tutorial. 

$> ./bin/mg5 
nig5> tutorial 

Then just follow the instructions on the screen and you will learn the basic command/usage 
of MadGraphS. This takes around 15-20 minutes. 

5.3 MC4BSM exercise 

In this section, we will study a full example on how to generate events for BSM theories. 
We will assume that you have your own UFO model [29] according to the FeynRules [1] 
tutorial. If you don't, you can download the associated model at the following address: 
http://feynrules.irmp.ucl.ac.be/attachment/wiki/WikiStart/MC4BSM_2012_UFO.tgzWe 
will split this example into three sections. First we will show how you can test the validity 
of the model, second, we will present how to create a valid set of parameters associated 
with the model, and finally we will present how you can generate BSM events, both with 
and without the associated decays. 

5.3.1 Importing and checking the model 

The simplest way to have access to a model in MG5 is to put it in the directory: MG5_DIR/models 
after that you can simply import it by typing 

ing5> import model MC4BSM_2012_UF0 

or 

mg5> import model MC4BSM_2012_UF0 --modelname 

The option -modelname tells MG5 to use the name of the particles defines in the UFO 
model, and not the usual MG5 conventions for the particles of the SM/MSSM. For this 
particular model, this changes only the name associated to r lepton (ta- and tt- respec- 
tively) . 

If you have developed your own model following the FeynRules tutorial from Sec- 
tion 3, this will be the first time that you are going to use this model. It is therefore crucial 
to start by checking the model. MadGraph5 performs some sanity checks the first time 
that you load the model, but those test are quite weak. We therefore suggest to test, three 
properties for on a series of processes. 

• The gauge invariance, by testing the Ward identities. 

• The Lorentz invariance. 

• The ALOHA [30] consistency, by evaluating the same square matrix element by 
different set of Helicity amplitudes. 
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For instance, we present how to check those properties for all the 2 — >^ 2 BSM particles 
productions: 

mg5> import model MC4BSM_2012_UF0 
mg5> define new = uv uv~ ev ev~ pi p2 
mg5> check p p > new new 

which results in the following output: 



Gauge results: 



Process 




matrix 




BRS 




ratio 




Result 


g 


g 


> uv 


uv~ 


7, 


.4113020914e- 


■01 


1, 


.0055761722e 


-31 


1. 


.3568144434e- 


-31 


Passed 


g 


u 


> uv 


pi 


3, 


.8373877873e- 


02 


6, 


. 1629758220e 


-33 


1, 


,6060341471e- 


-31 


Passed 


g 


u 


> uv 


p2 


2, 


.5726129500e- 


02 


4, 


.4176419953e 


-33 


1. 


,7171809678e 


-31 


Passed 


g 


u~ 


> uv~ pi 


2, 


.0117011717e- 


01 


8, 


.3456964257e 


-34 


4. 


. 1485766093e 


-33 


Passed 


g 


u~ 


■ > uv~ p2 


1, 


.9705216573e- 


01 


1, 


.8809915790e 


-32 


9. 


.5456529090e- 


-32 


Passed 



Summary: 5/5 passed, 0/5 failed 



Lorentz invariance results: 



Process 


Min element 




Max element 




Relative dif f . 




Result 


g 


g 


> uv uv~ 


4, 


.8743514998e- 


01 


4, 


. 87435 14998e 


-01 


1, 


, 1388417769e- 


16 


Passed 


g 


u 


> uv pi 


1, 


.3897148577e- 


•01 


1, 


.3897148577e 


-01 


1. 


. 1983282238e- 


15 


Passed 


g 


u 


> uv p2 


5, 


.7988729593e- 


02 


5, 


.7988729593e 


-02 


5. 


.9829676840e- 


16 


Passed 


g 


u~ 


> uv~ pi 


8, 


.8655875905e- 


03 


8, 


.8655875905e 


-03 


1. 


.7610238605e- 


15 


Passed 


g 


u~ 


> uv~ p2 


1, 


.0373598029e- 


01 


1, 


.0373598029e 


-01 


8. 


,0267932700e- 


16 


Passed 


u 


u 


> uv uv 


2, 


. 1536620557e+00 


2, 


. 1536620557e+00 


1. 


,8558171084e- 


15 


Passed 


u 


u~ 


> uv uv~ 


8, 


.4033626775e- 


01 


8, 


.4033626775e 


-01 


2. 


, 1138643036e- 


15 


Passed 


u 


u~ 


> pi pi 


7, 


.6001259005e- 


■03 


7, 


.6001259005e 


-03 


2. 


.9672410014e- 


15 


Passed 


u 


u~ 


> pi p2 


9, 


.5102790993e- 


05 


9, 


.5102790993e 


-05 


2. 


.7645774018e- 


14 


Passed 


u 


u~ 


> p2 p2 


4, 


.6427735529e- 


04 


4, 


.6427735529e 


-04 


8. 


.6404128750e- 


15 


Passed 


u 


u~ 


> ev ev~ 


2, 


.3919237366e- 


03 


2, 


.3919237366e 


-03 


2. 


,7196573767e- 


15 


Passed 


c 


c~ 


> uv uv~ 


8, 


.5752329113e- 


01 


8, 


.5752329113e 


-01 


2, 


,3304340127e- 


15 


Passed 


d 


d~ 


> uv uv~ 


8, 


.9275046680e- 


•01 


8, 


.9275046680e 


-01 


3. 


,2333557600e- 


15 


Passed 


d 


d~ 


> ev ev~ 


4, 


.8145327614e- 


04 


4, 


.8145327614e 


-04 


5. 


.6298410782e- 


16 


Passed 


s 


s~ 


> uv uv~ 


8, 


.8941849408e- 


01 


8, 


.8941849408e 


-01 


3. 


.6199458330e- 


15 


Passed 


u' 


u 


> uv~ uv~ 


2, 


.3800428911e+00 


2, 


. 380042891 le+00 


1. 


.8658874237e- 


15 


Passed 



Summary: 16/16 passed, 0/16 failed 

Not checked processes: c c~ > ev ev~, s s~ > ev ev' 



Process permutation results: 



Process 


Min element 




Max element 




Relative dif f . 




Result 


g 


g 


> uv uv~ 


4. 


,9184278197e- 


01 


4, 


.9184278197e 


-01 


2. 


.2572721717e- 


16 


Passed 


g 


u 


> uv pi 


4. 


,1859552985e- 


■02 


4, 


.1859552985e 


-02 


3. 


.3153215498e- 


16 


Passed 


g 


u 


> uv p2 


2. 


0129184319e- 


01 


2, 


.0129184319e 


-01 


1. 


, 1030978772e- 


15 


Passed 


g 


u~ 


> uv~ pi 


9. 


5566536137e- 


02 


9, 


.5566536137e 


-02 


5. 


,8086390357e- 


16 


Passed 


g 


u~ 


> uv~ p2 


3. 


,6165811126e- 


03 


3, 


.6165811126e 


-03 


3. 


.8372671248e- 


15 


Passed 
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u u > uv uv 2.1101603787e+00 2 . 1101603787e+00 2 . 1045282356e-15 Passed 

u u~ > uv uv~ 1.3549964258e+00 1 . 3549964258e+00 1 . 1470969258e-15 Passed 

u u~ > pi pi 4.9555140623e-03 4 . 9555140623e-03 1 . 4002369515e- 15 Passed 

u u~ > pi p2 2.0863569608e-02 2 . 0863569608e-02 6 . 6516842845e- 16 Passed 

u u~ > p2 p2 1.2914424155e-03 1 . 2914424155e-03 4 . 7013572521e- 15 Passed 

u u~ > ev ev~ 1 . 7823584674e-03 1 . 7823584674e-03 . OOOOOOOOOOe+00 Passed 

c c~ > uv uv~ 9.2608997797e-01 9 . 2608997797e-01 1 . 0789456164e-15 Passed 

d d~ > uv uv~ 8.3532448258e-01 8 . 3532448258e-01 1 . 3290919251e-16 Passed 

d d~ > ev ev~ 5 . 8126525280e-04 5 . 8126525280e-04 9 . 3262255680e-16 Passed 

u~ u~ > uv~ uv~ 2.1759606129e+00 2 . 1759606129e+00 2.0408880897e-16 Passed 
Summary: 15/15 passed, 0/15 failed 

More informations about these checks (hke the values of the random phase-space points) 
can be obtained via the commands: 

mg5> display checks 

The display commands is very useful in order to obtained information on the particles, 
couplings, interactions, ... For example^ 

mg5>display particles 

Current model contains 21 particles: 

ve/ve~ vm/vm~ vt/vt~ e-/e+ m-/m+ tt-/tt+ u/u~ c/c~ t/t~ d/d~ s/s~ 

b/b~ w+/w- uv/uv~ ev/ev~ a z g h pi p2 
mg5>display particles pi 

Particle pi has the following properties: 
{ 

'name' : 'pi' , 
'antiname' : 'pi ' , 
'spin' : 1, 
'color' : 1 , 
' charge ' : 0.00, 
'mass': 'MPel', 
'width': 'Wpel', 
'pdg_code': 9000006, 
'texname ' : 'pi ' , 
'antitexname ' : 'pi', 
'line' : 'dashed' , 
' pr opagat i ng ' : Tr ue , 
' is_part ' : True , 
'self_antipart' : True 

} 

mg5> 



^Note that the spin is written in the 25+1 convention. 
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5.3.2 Computation of the width and of the branching ratio 

The first requirement in order to have valid event generation is to have a vaUd set of 
input parameters. This includes having the correct widths and branching ratios (needed by 
Pythia). Currently, MadGraph5 is not able to compute these automatically and requires 
that you prescribe all decay channels. Note that this limitation will be resolve in a near 
future. 

So with the current version of MadGraph, you need to do the following in order to 
compute the width: 

nig5> import model MC4BSM_2012_UF0 
mg5> generate uv > u pi 
mg5> add process uv > u p2 
mg5> add process p2 > ev e+ 
mg5> add process p2 > ev~ e- 
mg5> add process ev > e- pi 
mg5> output 
mg5> launch 

The last command will ask you if you want to change the mass spectrum and/or the cuts. 
In principle all cuts should be set to zero for width computation. Then your browser will 
automatically open an html page with the results. ^'^ One of the results will be the new 
param_card.dat. Note that this script modified only the width of the particles present 
in the initial states of the generated processes. This means that the width of pi is not 
modified and therefore is still has it default value (1.0) which is not correct. You have to 
modify this card manually before using it for an event generation. 

5.3.3 Generating events 

In this section, we start the computation of the cross-sections and the generation of events 
for the proposed process of interest: 

pp — )• UU 

First we will generate this exact process, Pythia being in charge of the decays. Note that 
in this way, you lose the full spin-correlations. 

import model MODELNAME 
generate p p > uv uv~ 
output 
launch 

If you have install the PYTHIA-PGS/Delphes package, you will be asked if you want to run 
the package. Next a question will ask if you want to edit one of the cards. At this stage, 
you can enter the path of the card that you have created in the previous step and/or edit 
those manually. 

^"You can de-activated this behavior by modifying the configuration files present at the following path: 
./input / mg5 configuration.txt 
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In this case, the number of decay channels is quite hmited, it's therefore possible to 
prescribe all the steps of the decay chain. The syntax for the decay chains is the following: 
the production first, then the decays of the final states separate by a comma. To avoid 
ambiguity parenthesis should be present in the case of sub-decay. 

import model MC4BSM_2012_UF0 

generate p p > uv uv~ , uv > u pi, uv~ > u~ pi 
define 1 = e+ e- 
define Iv = ev ev~ 

add process p p > uv uv~, uv > u pi, \ 

(uv~ > u~ p2, (p2 > 1 iv, iv > 1 p2)) 
add process p p > uv uv~, uv~ > u~ pl,\ 

(uv > u p2, (p2 > 1 Iv, Iv > 1 p2)) 
add process p p > uv uv~, (uv > u p2, (p2 > 1 Iv, Iv > 1 p2)), \ 

(uv~ > u~ p2, (p2 > 1 Iv, Iv > 1 p2)) 

output 
launch 

The validity of this calculation must be considered, since specifying the decay sequence 
means the contribution of non-resonance Feynman diagram have been neglected. This is 
however valid since the interferences with those diagrams are negligible if the intermediate 
particles are on-shell. In order to ensure such a condition, we have associate to each of the 
decaying particles an additional cuts called BW_cut: 

{mmrt - mol < BWcut * r 

This cut can be specified in the run_card.dat. 

If you have installed the MadAnalysis package, you should have automatically dis- 
tributions created at parton level, after Pythia and at the reconstructed level -if runned-. 
Those one can be used to check the sanity of the productions but also to search the potential 
observables of this production. 

6 CalcHEP tutorial 

6.1 Installation of CalcHEP 
Download CalcHEP from: 

http : //theory. sinp .msu. ru/~pukhov/CALCHEP/calchep_3 .4 . . tgz 

Move it to the directory where you would like to do the tutorial. For example^^: 
mv ~/Downloads/calchep_3.4.0.tgz \ 

~/physics/pro j ects/MC4BSM-2012/Sandbox/calchep_3 .4.0. tgz 

Change to that directory. For example: 
^^If your OS already unzipped the file, change the endings appropriately. 
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cd ~/physics/projects/MC4BSM-2012/Sandbox 

Unpack the CalcHEP code^^: 
tar xvzf calchep_3.4.0.tgz 

Change to the codes dhectory: 
cd calchep_3.4.0 

Compile the code^'^: 
make 

Create a working directory for your calculations^^: 
. /mkUsrDir . . /ch_3 .4.0 

Change directory to the working directory: 
cd . ./ch_3.4.0 

If you have the XI 1 development libraries installed, you can start the symbolic session 
gui (otherwise skip this step): 

./calchep & 

Feel free to try some processes out and play. This will not be part of the tutorial. 

6.2 Install the MC4BSM 2012 model 

Download the CalcHEP model files from: 

http : //f eynrules . imp . ucl . ac . be/raw-attachment/wiki/WikiStart/MC4BSM_2012-CH . tgz 

Move this file to a convenient location. For example^^: 
mv ~/Downloads/MC4BSM_2012-CH.tgz \ 

~/physics/pro j ects/MC4BSM-2012/Sandbox/MC4BSM_2012-CH . tgz 

Change to this directory: 

cd ~/physics/projects/MC4BSM-2012/Sandbox 

'^^If your OS already unzipped, then you would do tar xvf calchep 3. 4.0. tar Or, if your OS already fully 
unpacked the file, you can skip this step. 

^^Note that you need the Xll development libraries installed on your system to use the graphical user 
interface (gui). For Mac, these used to come on the Xll disk. Now, I assume you get it through the app 
store. For Linux, your distribution likely has it in a repository. Google Xll development libraries and your 
distribution to see how to install it. However, note that it is not strictly necessary to have the gui working 
for this tutorial. 

^^mkUsrDir is planned to change to mkWorkDir. 

^^Again, change the ending as appropriate. 
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Unpack model files (if not already done): 
tar xzvf MC4BSM_2012-CH.tgz 

Install the model^^. We need to change the model number to one that is not already 
used. Since we just freshly installed CalcHEP, we know there are only three models 
included by default (the Standard Model with general CKM matrix, the Standard Model 
with diagonal CKM matrix and the Standard Model with effective Higgs vertices). The 
model number is given by the integer in the model file names. We will change it to 4. 

Change to the model file directory: 

cd MC4BSM_2012-CH 

Change the file names in such a way that the 1 becomes a 4. 

mv prtclsl .mdl prtcls4.mdl 
mv varsl.mdl vars4.nidl 
mv fund. mdl func4.mdl 
mv Igrngl.mdl lgrng4.mdl 
mv extlibl.mdl extlib4.mdl 

Now, copy the model files to the models subdirectory of the CalcHEP work directory: 
cp *4.mdl . ./ch_3.4.0/models/. 

Change directory back to the CalcHEP working directory: 
cd . ./ch_3.4.0 

If you have the XI 1 development libraries installed, you can start the symbolic session 
gui (otherwise skip this step): 

./calchep & 

You should see the tutorial model as the third model in the list. You can check your model 
for syntactical correctness, change it, calculate cross sections and generate distributions. 
Feel free to play. This will not be part of the tutorial. 

6.3 Generate events 

We, next, write a batch file to generate events. 

Open your favorite text editor (emacs, vi, gedit,...). For example: 

emacs & 

At the top of the file, we will choose the model. Add the lines: 

^®The model can be installed more easily using the gui, but I will describe the process that will work for 
everyone. 
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Model: 
Gauge : 



MC4BSM_2012 
Feynman 



The model name can be found at the top of the file models/prtclsS.mdl. 
We next specify the process we are interested in: 



Process : 


p,p->uv,uv~ 


Decay: 


uv->u,pl 


Decay: 


uv~->u~ ,pl 


Decay: 


uv->u,p2 


Decay: 


uv~->u~ ,p2 


Decay: 


p2->e- ,ev~ 


Decay: 


p2->e+,ev 


Decay: 


ev->e- ,pl 


Decay: 


ev~->e+,pl 


Alias: 


p=u,u~ ,d,d~ ,G 



We next specify the parton distribution function. We will use the built in CTEQ 6L 

set: 

pdfl: cteq61 (proton) 

pdf 2 : cteq61 (proton) 

Specify the momenta: 

pi: 4000 
p2 : 4000 

Specify the number of events to generate and the filename to use: 

Number of events (per run step) : 10000 
Filename : MC4BSM 

Specify your parallelization info: 

Parallelization method: local 
Max number of cpus : 2 
sleep time : 1 

I entered 2 for the max number of cpus because I have a dual core machine. You should 

change this as appropriate. 

Specify how many times Vegas should run: 

nSess_l: 5 

nCalls_l: 10000 

nSess_2: 5 

nCalls_2: 100000 
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The first five runs improve the grid. Then the statistics are cleared, the grid is frozen and 5 
more runs are done to calculate the cross section, widths and prepare the event generator. 

After you have all of these pieces entered, save the file. Call it "batch_file" in your 
work directory. 

Now, run the batch program by: 

. /calchep_batch batch_file 

The batch program will begin by stating where you can find information about the 
progress of your job. For example. 

Processing batch: 

Progress information can be found in the html directory. 
Simply open the following link in your browser: 

f ile : ///Users/neil/physics/pro j ects/MC4BSM-2012/Sandbox/ch_3 . 4 . 0/html/index . html 

You can also view textual progress reports in /Users/neil/physics/projects/MC4BSM-2012/Sand 

and the other .txt files in the html directory. 

Events will be stored in the Events directory. 

Note, the line that begins file:. You can copy and paste this line into your web 
browser's url bar. This shows the progress of your calculation. There is also a link on the 
left called Help which contains detailed instructions about how to write batch files. There 
are many options that we did not discuss in this tutorial. 

When the program finishes, you will have the following file: 

Events/MC4BSM-single . Ihe . gz 

which contains your events in LHE format. 
Enjoy! 

7 Pythia 8 tutorial 
7.1 Introduction 

This exercise corresponds to the Pythia 8 part of the more general MC tutorial given at 
the MC4BSM workshop at Cornell University 2012 [12]. It is for this reason focused on the 
particular task within this tutorial, however, should still serve as a good starting point to 
get familiar with the basics of how to use the PYTHIA 8 event generator. Much of these 
instructions were based on earlier PYTHIA 8 tutorials, which can be found on the homepage 
[31] and often include additional material than what is covered here. 

Within this first exercise it is not possible to describe the physics models used in the 
program; for this we refer to the online manual [32], PYTHIA 8.1 brief introduction [7], to 
the full Pythia 6.4 physics description [25], and to all the further references found in them. 

Finally, a good way to continue after the tutorial is often to chose a particular physics 
study your interest and then start to explore the different simulation and analysis aspects, 
using the different example programs together with the online manual, along the lines of 
the tutorial. 
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7.2 Installation and pre-workshop exercise 

Pythia 8 is, by today's standards, a small package. It is completely self-contained, and 
is therefore easy to install for standalone usage, e.g. if you want to have it on your own 
laptop, or if you want to explore physics or debug code without any danger of destructive 
interference between different libraries. 

Denoting a generic Pythia 8 version pythiaSlxx (at the time of writing xx = 62), 
here is how to install Pythia 8 on a Linux/Unix/MacOSX system as a standalone package. 

1. In a browser, go to 

http : / / www . thep .lu.se/ ~torb j orn/Pythia . html 

2. Download the (current) program package 

pythiaSlxx . tgz 

to a directory of your choice (e.g. by right-clicking on the link). 

3. In a terminal window, cd to where pythiaSlxx. tgz was downloaded, and type 

tar xvfz pythiaSlxx .tgz 
This will create a new (sub) directory pythiaSlxx where all the Pythia source files 
are now ready and unpacked. 

4. Move to this directory (cd pythiaSlxx) and do a make. This will take ~3 min- 
utes (computer-dependent). The Pythia 8 libraries are now compiled and ready for 
physics. 

5. For test runs, cd to the examples/ subdirectory. An Is reveals a list of programs, 
mainNN, with NN from 01 through 2S (and beyond). These example programs each 
illustrate an aspect of Pythia 8. For a list of what they do, see the "Sample Main 
Programs" page in the online manual (point 6 below). 

In order to test that the program installed successfully, most examples are suitable, 
however, for this tutorial (using version 8.162) a good choice would be mainll.cc, 
which reads in a test LHEF called ttbar.lhe. 
To execute this test program, do 

make main 11 

. /mainll . exe 

The output is now just written to the terminal, stdout, and if everything worked 
one should see (specific for mainll.cc) a text based histogram, showing the a charge 
particle multiplicity, at the end of the program output. To save the output to a file 
instead, do ./mainll. exe > mainll. out, after which you can study the test output 
at leisure by opening mainll. out. See Appendix A for a brief explanation of the 
event record. 

6. If you use a web browser to open the file 

pythiaSlxx/htmldoc/Welcome . html 
you will gain access to the online manual, where all available methods and parameters 
are described. Use the left-column index to navigate among the topics, which are then 
displayed in the larger right-hand field. 
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7.3 On-site exercise 



The task at hand for the Pythia 8 part of the tutorial is to read in the events, pro- 
vided in LHEF format from the earlier steps, and continue the event simulation following 
the hard scatter generation. Given that the LHEF includes both the information for ini- 
tialization and of the actual events, this is very straight forward. Copy first the LHEF 
events_with_ decay . Ihe .gz into the examples directory and unpack it using, 

gunzip event s_with_decay . Ihe . gz 
The LHEF should now be readable just like any ASCH file. 

The BSM model as well as the hard process is described in detail in the earlier steps of 
this tutorial, but as a short reminder the events consists of, pp — >• UU, production at the 
LHC with Ecm = 8 TeV. The U fermions have a mass of 500 GeV and the two following 
decay scenarios are possible: 

U^u^i] (7.1) 
U ^ u ^2 ^ u e E ^ u e e ^i. (7.2) 

The second step is to create our own main program. Open a new file mymain. cc in the 
examples subdirectory with a text editor, e.g. Emacs. Then type the following lines (here 
with explanatory comments added): 

// Headers and Namespaces. 

#include "Pythia. h" // Include Pythia headers, 

using namespace PythiaS; // Let PythiaS: : be implicit. 

int mainO { // Begin main program. 

// Set up generation. 

Pythia pythia; // Declare Pythia object 

pythia. readStringC "Beams : frame Type = 4"); // Beam info in LHEF. 
pythia. readStringC "Beams : LHEF = events_with_decay . Ihe") ; 
pythia. init ; // Initialize; incoming pp beams is default. 

// Generate event (s). 

for (int iEvent = 0; iEvent < 1; ++iEvent) { 
if (! pythia. next ) { 

if (pythia. info. atEndOf File 0) break; // Exit at enf of LHEF. 
continue; // Skip event in case of problem. 

} 

} 

pythia. Stat ; // Print run statistics, 
return 0; 

} 
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This will use one event from the LHEF and pass it through the remaining Pythia 8 
simulation. Go through the lines in the program and try to understand them by consulting 
the online manual [32]. 

Next you need to edit the Makefile (the one in the examples subdirectory) so it knows 
what to do with mymain.cc. The lines 

# Create an executable for one of the normal test programs 
mainOO mainOl main02 main03 . . . main09 mainlO mainlO \ 
and the three next enumerate the main programs that do not need any external libraries. 
Edit the last of these lines to include also mymain: 
mainSl . . . main40 mymain: \ 
Now it should work as before with the other examples: 
make mymain 

. /mymain. exe > mymain. out 

whereafter you can study mymain. out, especially the example of a complete event record 
(preceded by initialization information and by kinematical-variable listing for the same 
event). For this reason Appendix A contains a brief overview of the information stored in 
the event record. 

At this point you have in principle achieved your goal, the first event in the LHEF 
have been fully simulated by Pythia 8, using the default values of all settings. On the 
other hand, and despite only having one event, we only have access to the information of 
the individual particles in the event record, so interesting for technical validation but not 
much of real physics interest yet. 

7.4 A simple jet analysis 

We will now gradually expand the skeleton mymain program from above, in order to make 
a simple analysis, including jet reconstruction. The common standard, used in the experi- 
mental communities, is to produce fully simulated MC events in the HepMc format, either for 
analysis or further detector simulation. However, in order to keep the software infrastruc- 
ture at a minimum we will in this part only use the analysis functionality already available 
within Pythia 8. For further information regarding how to produce HepMc, the reader is 
referred to the online manual [32] together with other tutorials on the home page [31]. The 
final version of mymain.cc, with all steps included, is kept in Appendix B. Hence one can 
look there in case the instructions are unclear, however, that should of course be considered 
as the last resort, if getting completely stuck. 

The BSM events generated in this tutorial will always include jets, from the U decay, 
and missing $i particles. They will some times also contain electrons, depending on the 
U decay channel. For this reason we start with reconstructing jets from the final state 
particles in the event. For this we use the PYTHIA 8 Slow Jet program (found in the manual 
under Event Analysis), which is a simpler version of the commonly used Fast Jet jet finder. 
The SlowJet program supports both the kT, anti-kT, and Cambridge/ Aachen algorithms, 
where we will use the anti-kT which is a common choice within the LHC experiments. 
For further information about the jet reconstruction algorithms we recommend the online 
manual together with its references to the Fast Jet program. 



-34- 



Insert the following line before the event loop, in order to create a Slow Jet object, 

SlowJet: :SlowJet sJets(-l, 0.4, 50., 5., 1, 2); 
The first argument specifies the anti-kT algorithm, the second is the related R measure 
which roughly corresponds to the radius of the jet cone in (y, (?!)). The next two corresponds 
to, minimum pT and maximum \ri\, acceptance cuts usually implied by the experiment of 
interest. The last two arguments specifies respectively to use all final state particles and 
to use their actual mass in the reconstruction. The selection of which particles to include 
is normally done more carefully, however, in a trade-off between illustrating the physics 
without spending too much time on technical details and infrastructure, we chose a rather 
unsophisticated approach here. Just after this line, you should also create the following two 
histograms. 

Hist nJ( "Nr Jets", 10, -0.5, 9.5); 

Hist jlpT( "Leading Jet pT", 100, 0., 1000.); 
which corresponds to our final analysis results. 

The next step is to reconstruct jets from the final state particles in each event according 
to the jet algorithm we just defined. After that we also want to fill our histograms, once 
per event, and at the end of the program print them to our output. To do this, include 
the following lines inside the event loop, but at the end so that pythia.nextO has been 
executed, 

// Jet analysis. 

if ( s Jet s . analyze (pythia. event) ) { 
nJ.filKsJets.sizeJetO) ; 
jlpT.filKsJets.pKO)) ; 

} 

The pythia member . event corresponds to the event record and you can try to find out the 

meaning of these lines by using the online manual. To write the histograms to the output 
include the following line at the very end, just before the return statement, 
std: :cout « nJ « jlpT; 

Now increase the number of events in the event loop, e.g. to 1000, compile and run the 
program. At the end of the output one should now see the number of jets and leading jet 

Pt distributions, in traditional text based histograms. 

During the run you may receive problem messages. These come in three kinds: 

• a warning is a minor problem that is automatically fixed by the program, at least 
approximately; 

• an error is a bigger problem, that is normally still automatically fixed by the program, 
by backing up and trying again; 

• an abort is such a major problem that the current event could not be completed; in 
such a rare case pythia.nextO is false and the event should be skipped. 

Thus the user need only be on the lookout for aborts. During event generation, a problem 
message is printed only the first time it occurs. The above-mentioned pythia. statO will 
then tell you how many times each problem was encountered over the entire run. 
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What we did so far is clearly not what we want to do, since first we include the $i 
which should be invisible and second the electrons produced in the decay chains, potentially 
also with high pT, are included. 

For the purpose of this tutorial we use a rather ugly trick in order to eliminate particles 
with respect to the jet reconstruction. We simply loop through the event record and when we 
find a $i (id == 9000006) we set its status to a negative value, since the jet reconstruction 
only considers final state particles, i.e. with positive status. In order to do this, insert the 
following lines, just before the jet analysis part, 

for (int iPart = 0; iPart < pythia. event . size () ; ++iPart) { 
if (pythia. event [iPart] .idAbsO == 9000006) { 
int Stat = pythia. event [iPart] . status ; 
if (stat > 0) pythia. event [iPart] . status(-stat) ; 

} 

} 

These lines hence illustrate both how to loop over the particles in the event record and how 
to access their properties, in this case the id and status codes. Now compile and run the 
program again. The mean number of jets should now have been reduced by approximately 
two (from 4.8 to 3.1), due to the two <I>i in the events. The leading px jet should now 
dominantly originate from the quarks in the U decay and the mean pT therefore is around 
half of the U mass, 500 GeV. 

We now leave it as an exercise to also eliminate any final state electrons in the event 
and see how the distributions are affected. Another easy thing to investigate would be to 
see how the number of reconstructed jets varies with the minimum pr requirement specified 
for the jet finder. 

7.5 The event record 

The event record is set up to store every step in the evolution from an initial low-multiplicity 
partonic process to a final high-multiplicity hadronic state, in the order that new particles 
are generated. The record is a vector of particles, that expands to fit the needs of the current 
event (plus some additional pieces of information not discussed here). Thus event [i] is 
the i'th particle of the current event, and you may study its properties by using various 
event [i] .methodO possibilities. 

The event .listO listing provides the main properties of each particles, by column: 

• no, the index number of the particle (i above); 

• id, the PDG particle identity code (method idO); 

• name, a plaintext rendering of the particle name (method nameO), within brackets 
for initial or intermediate particles and without for final-state ones; 

• status, the reason why a new particle was added to the event record (method 
status 0); 

• mothers and daughters, documentation on the event history (methods motherlO, 
mother2(), daughterlO and daughter2()); 
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• colours, the colour flow of the process (methods col() and acolO); 

• p_x, p_y, p_z and e, the components of the momentum four-vector {px,Py,Pz, E), in 
units of GeV with c = 1 (methods px(), pyO, pz() and e()); 

• m, the mass, in units as above (method mO). 

For a complete description of these and other particle properties (such as production and 
decay vertices, rapidity, p±, etc), open the program's online documentation in a browser 
(see Section 7.2, point 6, above), scroll down to the "Study Output" section, and follow the 
"Particle Properties" link in the left-hand-side menu. For brief summaries on the less trivial 
of the ones above, read on. 

7.5.1 Identity codes 

A complete specification of the PDG codes is found in the Review of Particle Physics [33]. 
An online listing is available from 

http : / /pdg . Ibl . gov/2008/nicdata/nic_particle_id_contents . html 
A short summary of the most common id codes would be 



1 


d 


11 


e 


21 


g 


211 




111 




213 




2112 


n 


2 


u 


12 




22 


7 


311 


KO 


221 


V 


313 


K*o 


2212 


P 


3 


s 


13 




23 


ZO 


321 


K+ 


331 


v' 


323 


K*+ 


3122 


AO 
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c 


14 




24 


W+ 


411 


D+ 


130 


K 


113 




3112 


s- 


5 


b 


15 




25 


HO 


421 


DO 


310 


Kg 


223 




3212 




6 


t 


16 








431 








333 





3222 





Antiparticles to the above, where existing as separate entities, are given with a negative 
sign. 

Note that simple meson and baryon codes are constructed from the constituent (anti) quark 
codes, with a final spin-state-counting digit 2s -|- 1 (K^ and Kg being exceptions), and with 
a set of further rules to make the codes unambiguous. 

7.5.2 Status codes 

When a new particle is added to the event record, it is assigned a positive status code that 
describes why it has been added, as follows: 



code 


range 


explanation 


11 


- 19 


beam particles 


21 


- 29 


particles of the hardest subprocess 


31 


- 39 


particles of subsequent subprocesses in multiparton interactions 


41 


- 49 


particles produced by initial-state-showers 


51 


- 59 


particles produced by final-state-showers 


61 


- 69 


particles produced by beam-remnant treatment 


71 


- 79 


partons in preparation of hadronization process 


81 


- 89 


primary hadrons produced by hadronization process 


91 


- 99 


particles produced in decay process, or by Bose-Einstein effects 



Whenever a particle is allowed to branch or decay further its status code is negated (but it 
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is never removed from the event record), such that only particles in the final state remain 
with positive codes. The isFinalO method returns true/false for positive/negative 
status codes. 

7.5.3 History information 

The two mother and two daughter indices of each particle provide information on the history 
relationship between the different entries in the event record. The detailed rules depend on 
the particular physics step being described, as defined by the status code. As an example, 
in a 2 — )• 2 process ab — )■ cd, the locations of a and b would set the mothers of c and d, with 
the reverse relationship for daughters. When the two mother or daughter indices are not 
consecutive they define a range between the first and last entry, such as a string system 
consisting of several partons fragment into several hadrons. 

There are also several special cases. One such is when "the same" particle appears as 
a second copy, e.g. because its momentum has been shifted by it taking a recoil in the 
dipole picture of parton showers. Then the original has both daughter indices pointing to 
the same particle, which in its turn has both mother pointers referring back to the original. 
Another special case is the description of ISR by backwards evolution, where the mother 
is constructed at a later stage than the daughter, and therefore appears below in the event 
listing. 

If you get confused by the different special-case storage options, the two pythia . event . motherList ( i) 
and pythia. event . daughterList (i) methods are able to return a vector of all mother 
or daughter indices of particle i. 

7.5.4 Color flow information 

The colour flow information is based on the Les Houches Accord convention [34]. In it, the 
number of colours is assumed infinite, so that each new colour line can be assigned a new 
separate colour. These colours are given consecutive labels: 101, 102, 103, .... A gluon 
has both a colour and an anticolour label, an (anti)quark only (anti)colour. 

While colours are traced consistently through hard processes and parton showers, the 
subsequent beam-remnant-handling step often involves a drastic change of colour labels. 
Firstly, previously unrelated colours and anticolours taken from the beams may at this stage 
be associated with each other, and be relabeled accordingly. Secondly, it appears that the 
close space-time overlap of many colour fields leads to reconnections, i.e. a swapping of 
colour labels, that tends to reduce the total length of field lines. 

7.6 The jet analysis program 

This is the final mymain. cc at the end of the tutorial, 

#include "Pythia. h" // Include Pythia headers. 

using namespace PythiaS; // Let PythiaS:: be implicit. 

int mainO { // Begin main program. 
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// Set up generation. 

Pythia pythia; // Declare Pythia object 

pythia.readStringC "Beams :frameType = 4"); // Beam info in LHEF. 
pythia. readStringC "Beams :LHEF = events_with_decay .Ihe") ; 
pythia. init ; // Initialize; incoming pp beams is default. 

SlowJet: :SlowJet sJets(-l, 0.4, 50., 5., 1, 2); 

Hist nJ( "Nr Jets", 10, -0.5, 9.5); 

Hist jlpT( "Leading Jet pT" , 100, 0., 1000.); 

// Generate event(s). 

for (int iEvent = 0; iEvent < 1000; ++iEvent) { 
if (! pythia. next ) { 

if (pythia. info.atEndOf File 0) break; // Exit at enf of LHEF. 
continue; // Skip event in case of problem. 

} 

// Loop through event record. 

for (int iPart = 0; iPart < pythia. event . size () ; ++iPart) { 
if (pythia. event [iPart] .idAbsO == 9000006) { 
int Stat = pythia. event [iPart] . status ; 
if (stat > 0) pythia. event [iPart] . status (-Stat) ; 

} 

if (pythia. event [iPart] .idAbsO == 11) { 
int Stat = pythia. event [iPart] .status ; 
if (stat > 0) pythia. event [iPart] . status (-Stat) ; 

} 

} 

// Jet analysis. 

if (s Jets . analyze (pythia. event) ) { 
nJ.filKsJets.sizeJetO) ; 
jlpT.filKsJets.pKO)) ; 

} 

} 

pythia. stat ; // Print run statistics. 
std::cout « nJ « jlpT; 
return 0; 
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8 HERWIG++ tutorial 

8.1 Introduction 

HerwigH — h [35, 36], is a general-purpose Monte Carlo event generator for the simulation 
of hard lepton-lepton, lepton-hadron and hadron-hadron collisions. 
Some of the main features, in random order, are: 

1 . Many matrix elements natively and capability of reading in Les Houches-accord event 
files. 

2. A description of the underlying event via multiple parton interactions. 

3. Initial- and final-state QCD parton showering. 

4. Next-to-leading order: Several POWHEG [37, 38] processes and compatibility with 
MC@NLO [39]. 

5. Cluster hadronization. 

6. Individually modelled hadron and tau decays. 

7. QED radiation. 

8. Analysis via RiVET [40], native analysis handlers, or output to HepMC events. 

This tutorial aims to guide you through the steps of going from the HERWIG++ source 
distribution through to some basic analysis for a specific model, which is described in 
Section 2. The model events are generated in MadGraph [22], and will be 'fed' into the 
event generator via the Les Houches interface^^. 

Note that Wiki page for this tutorial also exists at: 
http : / /herwig. hepf orge . org/trac/wiki/MC4BSM. 

8.2 Installation and setup 
8.2.1 Required libraries 

Throughout this tutorial, it will be assumed that you have created a working directory in 
your home directory, which we will be referring to as $MCBSM from now on. If you wish, you 
may define $MCBSM (in bash shell): 

export MCBSM=/target/directory 

In c-shell you need to use: 

setenv MCBSM "/target/directory" 

All the source tarballs will be assumed to be placed in $MCBSM/source. 

If the libraries do not already exist on your system, you will need to compile and install 
the following: 

^''in the future, Herwig++ will include automated creation of the model internally via a FeynRules 
UFO interface [29]. 
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• GNU Scientific Library (GSL), 

http: //mirror . switch. ch/f tp/mirror/gnu/gsl/gsl- 1 . 15 . tar . gz. 

Recommended installation in $MCBSM/gsl. 

• LHAPDF, the Les Houches Accord PDF Interface, 

http: //www. hepf orge . org/archive/lhapdf /lhapdf -5 . 8 . 7 . tar . gz. 
Recommended installation in $MCBSM/ lhapdf . 

• The PDF used in this tutorial, cteqGll, 

http: //www. hepf orge . org/archive/lhapdf /pdfsets/current/cteq611 . LHpdf . 
To be placed in $MCBSM/lhapdf /share/lhapdf /PDFsets (the directory may need to 
be created). 

Important: Make sure that any program binaries and libraries you install appear in your 
$PATH and $LD_LIBRARY_PATH$ ($DYLD_LIBRARY_PATH in MacOS X) respectively, otherwise 
the binaries or libraries will not be found by the system. 

Note that when compiling, if you have a multi-core machine, you may use make -jX, 
where X is the number of cores, to speed up the process. 

GSL Compilation and installation in the conventional way: 

tar xvzf gsl-1 . 15.tar .gz; cd gsl-1.15 
./configure --pref ix=$MCBSM/gsl 
make 

If everything is successful (i.e. there are no errors): 
make install 

LHAPDF Equivalent to GSL with the necessary change of target directory through: 
— pref ix=$MCBSM/lhapdf 

Once installed, make sure cteqGll. LHpdf is placed in $MCBSM/lhapdf /share/lhapdf /PDFsets. 
8.2.2 Installing HERWIG++ and ThePEG 

HERWIG++ is based on ThePEG, (Toolkit for High Energy Physics Event Generation), 
which is essentially a library that contains the necessary building blocks for constructing 
an event generator. 

First we need to compile and install ThePEG. Get the latest version at: 
http : / /www. hepf orge . org/archive/thepeg/ThePEG- 1.7. 3 . tar .bz2. 
To compile and install: 

tar xvfj ThePEG- 1.7. 3. tar. bz2; cd ThePEG- 1.7. 3 
./configure --pref ix=$MCBSM/thepeg --with-gsl=$MCBSM/gsl 
— with-LHAPDF=$MCBSM/lhapdf 
make; make install 



- 41 - 



For instructions on how to do this on the MacOS X operating system, consult: 
http : //herwig. hepf orge . org/trac/wiki/HwOsX 

Once ThePEG is installed, you are ready to install HerwigH — h, linking to all the 
necessary libraries. Get the latest stable release (currently version 2.5.2) at: 
http : //www. hepf orge . org/archive/herwig/Herwig++-2 . 5 . 2 . tar .bz2. 
To compile and install: 

tar xvfj Herwig++-2.5.2.tar.bz2; cd Herwig++-2 . 5 . 2 
./configure --pref ix=$MCBSM/hpp --with-thepeg=$MCBSM/thepeg 
--with-gsl=$MCBSM/gsl 
make; make install 

8.2.3 HERWIG++ preliminaries 

To run HerwigH — h one needs to set up a working directory. In this tutorial, this will be 
taken to be $MCBSM/hppwork. Create this directory and copy the following: 

cp $MCBSM/source/Herwig++-2.5.2/src/*.in $MCBSM/hppwork 

cp $MCBSM/source/Herwig++-2.5.2/src/*.model $MCBSM/hppwork 

cp $MCBSM/source/Herwig++-2.5.2/src/defaults/*. in $MCBSM/hppwork 

This copies all the necessary input files, including the ones containing the default settings, 
as well as different .model files (such as the MSSM, ADD, and so on) for future use. Once 
this is done, change into the $MCBSM/hppwork directory and execute: 

Herwig++ init 

If the program has compiled correctly, the binaries and libraries are in your $PATH and 
$LD_LIBRARY_PATH respectively, and all the necessary . in files have been copied properly, 
there should be no errors. 

Most of the copied .in files contain inputs for the default settings (e.g. PDF. in in- 
cludes defaults for PDF handling. Cuts . in includes the default cuts and so on). It is not 
recommended in general to change the settings for running in those files directly. Instead, 
one should use a customized input file for the individual runs. The simplest example to 
look at is the file LHC . in. In there, you need to choose the desired energy at which you 
wish to run the event generator. The default in this input file is LHC at 7 TeV: 

######################## 
## sqrt(s) = 7000 GeV ## 
######################## 

set LHCGenerator :EventHandler :LuminosityFunction:Energy 7000.0 
set /Herwig/UnderlyingEvent/KtCut :MinKT 3.06 
set /Herwig/UnderlyingEvent/UECuts:MHatMin 6.12 

^*The subdirectory source will not be required if you have not put the source files in there! 
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You may wish to examine the construction of the LHCGenerator in the file HerwigDef aults . in, 
which is of type EventHandler.^^ . It is recommended not to make any changes in 
HerwigDef aults . in, instead, these should be made in LHC.in, as also mentioned previ- 
ously. A useful thing to be able to do is change the type of the hard process to be generated. 
The default in LHC.in is Drell-Yan Z/photon: 

# Drell-Yan Z/gamma 

insert SimpleQCD:MatrixEleinents [0] MEqq2gZ2ff 

Let's comment this out (by adding a # in front) and instead uncomment: 

# top-antitop production 

insert SimpleQCD:MatrixEleinents [0] MEHeavyQuark 
and execute: 
Herwig++ read LHC.in 

If there are no errors, the command should return no text output. It should create a file 
LHC.run, which will contain internal instructions to the event generator, which are not 
human-readable. To make a short run, say, of 100 events of top-antitop events, including 
initial/final-state radiation, hadronization, and the underlying event, simply type: 

Herwig++ run LHC.run -NlOO 

The run will initialize and commence generating events. This point is a good place to 
stop the pre-workshop part of the tutorial. 

Running the file with the above command will not produce any significant output apart 
from the file LHC.out, which will contain the cross section of the chosen process (you may 
even choose more than one process). A log file may obtained by adding -dl to the above: 

Herwig++ run LHC.run -NlOO -dl 

This will print the first 10 events of the run (determined by setLHCGenerator : PrintEvent 10 

in the LHC . in file) in LHC . log. Have a look at it and see if you can recognize the different 

steps of the event generation, and follow through the decay products of the top and antitop. 

For a full explanation have a look at the Wiki page: 

http : //herwig. hepf orge . org/trac/wiki/EventRecordFormat. 

You may also wish to enable the built-in analysis for top-antitop: 

# analysis of top-antitop events 

insert LHCGenerator: AnalysisHandlers /Herwig/ Analysis/TTbar 

Note that every time a change is made in the . in file, one needs to invoke the 'read' 
command. 

After running the event generator again, the analysis will produce a file LHC-TTbar . top. 
This is made to be read by TopDrawer, a rather archaic, but useful, plotting package. 

^®The EventHandler class is a general class of ThePEG, responsible for the generation of events. 



-43- 



Unfortunately it is difficult to compile it from source, and pre-compiled binaries may or 

may not work on your system. You may try the binaries subdirectory of: 

http : / /www . rcnp . osaka- u . ac . jp/~okamura/ ftp/ pub/ index . php?dir=topdrawer/ 

If the executable works, you may try to create a postscript file by running: 

td -d ps LHC-TTbar, 

and using Ghost View (gv) to view it. 

To run through the rest of the tutorial, you will need to obtain the following tarball: 
http : //www. itp .uzh. ch/~andreasp/mc4bsm. tar . gz 

Extracting this in $MCBSM should create a subdirectory mc4bsm. The details of the contents 
of $MCBSM will be discussed in Section 8.3.1. 

8.3 The on-site tutorial 

This tutorial will study the toy model described in Section 2. We assume the existence 
of MadGraph parton-level (but with the heavy new particles already decayed) event files 
(see Section 5). 

8.3.1 Running HERWIG++ with the toy model 

Extracting inc4bsin.tar.gz in $MCBSM should create a subdirectory mc4bsin. This in turn 
contains the following subdirectories: 

• analysis: contains the sources of a basic HerwigH — h analysis. 

• hppinput: contains the input files necessary to run the basic tutorial events. 

• madgraph_events: contains the MadGraph events which include decays of the new 
heavy particles down to 'stable' particles. 

• results: contains 'solution' plots to the tutorials. 

We first need to compile the customized analysis handler. To do that, copy the 
MC4BSMAnalysis.tar.gz to the Herwig++-2.5.2/Contrib subdirectory, and extract: 

tar xvzf MC4BSMAnalysis.tar.gz 

This should overwrite the default Herwig++-2 . 5 . 2/Contrib/Makef ile . am and put the 
analysis source code in the MC4BSMAnalysis. In the Herwig++-2 . 5 . 2/Contrib directory, 
type: 

sh inake_inakef iles . sh 

This forces the recreation of all the Makefiles from input, and hence also creates a Makefile 
for our analysis. Change to the Herwig++-2 . 5 . 2/Contrib/MC4BSMAnalysis directory and 
type make. This should compile the library MC4BSMAnalysis . so. Copy this into the Her- 
WIGH — h installation directory library: 

cp MC4BSMAnalysis. so $MCBSM/hpp/lib/Herwig++ 
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In your HerwigH — h working directory ($MCBSM/hppwork), copy the contents of hppinput 
and niadgraph_ events. These should be the input file MC4BSM.in and the model file 
MC4BSM . model from the prior directory and event s_with_decays . Ihe from the latter. Have 
a look and try to understand at what's in each of the MC4BSM. * files: The MC4BSM. in file is 
equivalent to the LHC . in input file, with the appropriate settings to read in a Les Houches- 
accord event file. The underlying event settings have been chosen for an anticipated 8 TeV 
run. There's also a change of PDF to the one used by the input events (cteqGU), using the 
LHAPDF library. Even though the input event file provides all the necessary information for 
the hard process, HerwigH — h needs to know about the existence of the new particles: that 
is why the second file, MC4BSM. model, is necessary, and is read at the top of MC4BSM.in. 
It's easy to construct one manually:^'^ for example, the following defines particle uv with 
PDG id 9000008, mass 500 GeV and width 1 GeV. 

create ThePEG: :ParticleData uv 

setup uv 9000008 uv 500.0 1.0 10.0 1 . 973269631e-13 2 3 2 

The rest of the numbers in order are: maximum number of widths away from central value 
of mass allowed, the lifetime (calculated from the width), the electromagnetic charge (units 
of 1/3), the SU{3)c charge (3 implies triplet, 1 implies singlet and so on), the spin 2s + 1 
and a boolean for whether the particle is stable or not. 
Read in and run the input file in the usual way: 

Herwig++ read MC4BSM.in 
Herwig++ run MC4BSM.run -N9000 

The analysis should create a subdirectory mcplots containing some plots, in TopDrawer 
format. You may view them using TopDrawer, or copy the data in your favourite plotting 
program. The plot files are the pseudorapidities and transverse momenta of the first 6 
highest-pr jets in the event (where in this case we are also including the invisibles), the 
number of jets (Njets) and the total invariant mass in the event (including the invisibles), 
calculated by: 

here E and P are the sums of the all the individual particle energies and 3-momenta 
respectively, within a certain range of pseudorapidity. The pseudorapidity cut for the jets, 
as well as other options for the analysis can be found at the relevant section in the input 
file: 

# Choose jet clustering algorithm 

set /Herwig/Analysis/MC4BSMAnalysis : JetAlgorithm 2 

# Jet ET cut to apply in jet clustering 

set /Herwig/Analysis/MC4BSMAnalysis :ETClus 20*GeV 

# Cone size used in clustering in merging, 
set /Herwig/Analysis/MC4BSMAnalysis :RClus 0.4 

^"in a future release, it will be possible to construct it using the FeynRules interface. 
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# Max I eta I for jets in clustering in merging, 
set /Herwig/Analysis/MC4BSMAnalysis :EtaClusMax 5 

The options ETClus, RClus and EtaClusMax correspond to the minimum pT of the jets, the 
clustering parameter for the algorithm and the maximum pseudorapidity for the jets. 

Try and see what happens when you switch off the hadronization and underlying event, 
by uncommenting the relevant pieces: 

#set LesHouchesHandler :HadronizationHandler NULL 
#set /Herwig/Shower/ShowerHandler :MPIHandler NULL 

You may wish to edit the analysis code and see if you can add an extra plot. This can be 
done by editing MC4BSMAnalysis . cc and MC4BSMAnalysis . h in the Contrib/MC4BSMAnalysis 
subdirectory. You will need to add a pointer for the new histogram in the . h file and then 
add the analysis code in the . cc file to declare the limits of the histogram, fill it in, and 
write it out. Take the existing histograms as an example. A good one to try is the pT of all 
leptons, for which there already exists commented-out code at the right places in the code 
for guidance: search for "EXAMPLE NEW HISTOGRAM:" in the code. There should be three 
points in the . cc file and one in the . h file. 

8.3.2 Further reading 

For specifics of the HerwigH — h event generator, the manual and release notes for version 
2.5 are useful [35, 36]. For further, up-to-date information, you may consult the Wiki pages, 
at: 

http : / /herwig. hepf orge . org/trac/. 

For details on the specific inner workings of HerwigH — h and ThePEG, consult the 
Doxy gen output at: 

http : / /herwig . hepf orge . org/ doxygen/ index . html 
http : //thepeg. hepf orge . org/doxygen/ index. html 

Although it may look a bit difficult to navigate through, it is a good starting point to find 
out how things are interconnected inside the event generator. 

9 Sherpa tutorial 
9.1 Requirements 

You will need a PC running Linux or a Mac to follow this exercise. Many of the Macs the 
Sherpa team came across required some special tweaks to guarantee smooth installation. 
Linux should therefore be the preferred operating system for this exercise. 

You should have the GNU autotools installed. They will also come in handy for other 
applications. Try to locate libtoolize, automake and autoconf on your system. If they are 
missing, and there is no package manager that provides easy installation, you can download 
the sources from 

• http:/ /ftp. gnu.org/gnu/autoconf/autoconf-2.65.tar.gz 



-46- 



• http:/ /ftp. gnu. org/gnu/automake/automake-1. 11.1. tar. gz 

• http:/ /ftp.gnu.org/gnu/libtool/libtool-2.2.6b.tar.gz 

Follow the installation instructions provided in the tarballs. 

Although not mandatory for this exercise, it will be useful to have metapost installed. 
It will be used to plot the Feynman graphs generated by Sherpa. 

To visualize the analysis results, you will need gnuplot or a similar plotting tool. The 
instructions given here assume that you have gnuplot installed. 

9.2 Installation 

There are two methods to install Sherpa on your system 

9.2.1 ... via package download 

The Sherpa sources for this exercise are obtained from either of these sites 

• http: / /sherpa. hepforge.org/trac/wiki/SherpaDownloads/Sherpa-1. 4.0 

• http:/ /www. slac.stanford.edu/~shoeche/mc4bsm/Sherpa-l. 4.0 

After the package download go to the download directory and follow these steps 

tar -xzf SHERPA-MC- 1.4.0. tar. gz 
cd SHERPA-MC- 1.4.0/ 
./configure --enable -analysis 
make -j2 install 

If successful, you should now find a Sherpa executable in the bin/ subdirectory. On a Mac 
with 64 bit system, use FC='gf ortran -in64' . /configure --enable-analysis instead of 
the above configure command. 

9.2.2 ... via svn checkout 

If you have a subversion client installed on your system, you can check out Sherpa as 
follows. This will be an easier option if you try to keep your version up-to-date with bugfixes. 
Note that you must have the GNU autotools installed to proceed with this method. Follow 
these steps 

svn CO http://sherpa.hepforge.org/svn/branches/rel-l-4-0 
cd rel-1-4-0/ 
autoreconf -i 

./configure --enable-analysis 
make -j2 install 

If successful, you should now find a Sherpa executable in the bin/ subdirectory. On a Mac 
with 64 bit system, use FC='gf ortran -m64' . /configure --enable-analysis instead of 
the above configure command. 
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9.2.3 Dealing with problems 

If there are problems during the installation process, please send an email to info@sherpa- 
mc.de or list the issue on our bug tracker. Some advice can also be found in Section 2 - "Get- 
ting started" - of the online manual, which is available at http: / /sherpa.hepforge.org/doc/SHERPA- 
MC-1.4.0.html 

9.3 Testing your Sherpa installation 

To guarantee a successful installation of Sherpa, we will try to simulate the production of 
T^^-boson plus jets final states at the Tevatron using ME+PS merging with up to two jets 
only (for simplicity). Follow these steps 

cd Examples/Tevatron_WJets/ 
. . / . . /bin/Sherpa 

The example should take about 15 minutes to run. 
You should see the following output 

Welcome to Sherpa, <useriiame>. Initialization of framework underway. 
Run_Parameter: : InitO : Setting memory limit to <limit> GB. 



Event generation run with SHERPA started 



... lots more following ... 

+ 

I 

I Total XS is 2338.05 pb +- ( 22.9361 pb = 0.98 7, ) 
I 

+ 

... some more following ... 

Time: 17m 51s on Thu Mar 8 12:31:59 2012 
(User: 17m 49s, System: 2s, Children User: Os, Children System: Os) 

If this output shows up in your terminal, you have successfully installed Sherpa. 

Congratulations! 
9.3.1 Dealing with problems 

If Sherpa refuses to run, please send an email to info@sherpa-mc.de or list the issue on 
our bug tracker. 
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9.4 Getting started 

The setup is available from 

• http:/ /www. slac.stanford.edu/~shoeche/mc4bsm/2012/ 
Download the tarball mc4bsni_sherpa.tar.gz and unpack the setup files using 
tar -xzf inc4bsin_sherpa.tar.gz 

Note that we use a given set of FeynRules output files, which are stored in the sub- 
directory FeynRules_Output. You can overwrite these files with your own. However, if you 
alter the masses of U, E, (pi or (j)2, the particle widths have to be recomputed. Section 9.8 
explains how to do this. You can work through Sec. 9.8 first, but it is recommended to 
follow Sees. 9.5-9.8 in order, as this will help to understand the input structure of Sherpa. 

For the remainder of this exercise we will assume that Sherpa was installed into a 
directory called <sherpa_pref ix>/. You should be able to locate the executeable in the 
subdirectory <sherpa_pref ix>/bin/, cf. the pre-workshop instructions. 

9.5 Understanding the input structure 

Change to the directory Intro and open the input file Run.dat. The file consists of various 
sections, which are marked as (model), (beam), (process), etc. 

The (model) section contains an instruction to use the SHERPA-internal FeynRules 
interface and to flag the particles E, (pi and (p2 as unstable. 

The (beam) section is used to set up the collider type {pp) and its cms energy (8 TeV). 

The (processes) section defines the reaction of interest. Here we simulate the process 
uu — )• UU. We also instruct Sherpa to write out Latex files depicting the Feynman graphs. 

The (me) section is used to set a scale at which the strong coupling is to be evaluated. 

Run the simple example using 

<sherpa_pref ix>/bin/Sherpa 

Sherpa will stop with the message 

New libraries created. Please compile. 

followed by some citation info. Now you need to compile the process-specific source code 
generated by Sherpa using 

. /makelibs 

After the compilation has finished, run Sherpa again 

<sherpa_pref ix>/bin/Sherpa 

The program will now dynamically link the libraries which have just been created and 

compute a cross section for the process uu — )• UU. 

Have a look at the Feynman graphs that contribute to this process 

. /plot_graphs Process/P2_2/ 
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9.6 Simulating parton-level events 



Change back to the original directory and then go to ToyModel_PartonLevel. 
Have a look at the Run . dat input file. 

In the (run) section we disable the hadronzation module by using FRAGMENTATION Off, 
the parton shower by using SHOWER_GENERATOR None and the YFS soft photon generator 
by using ME_QED Off. 

In the (processes) section, we instruct Sherpa to generate the following reactions: 

• pp ^ U[-^ (piu] (piu] 

• pp ^ U[-^ 02 [—^ 01 e^e^] u] )• 

• pp ^ U[-^ )• 02 [—J" 01 e+e^] u\ 

• pp ^ [/[— > 02[— > 01 e+e"] u] 02 [-^ 01 e+e"] u] 
Run this setup using 

<sherpa_pref ix>/bin/Sherpa 

Again, Sherpa will stop with the message 

New libraries created. Please compile. 

followed by some citation info. You need to compile the process-specific source code gen- 
erated by Sherpa using ./makelibs. After the compilation has finished, run Sherpa 
again 

<sherpa_pref ix>/bin/Sherpa 

The cross sections will be computed and Sherpa will generate 10000 events. Near the end 
of the output you should see the following message 

+ 

I 

I Total XS is 0.444465 pb +- ( 0.00429378 pb = 0.96 7, ) 
I 

+ 

where the precise value of the cross section depends on the mass parameters you have chosen 
for the [/, E and 0^/2 fields. It is important that you use this cross section, when computing 
event rates, rather than any of the cross sections quoted during the integration step. A 
detailed explanation why can be found in section three of the Sherpa online manual, 
http:/ /sherpa. hepforge.org/doc/SHERPA-MC-1. 4. 0.html^Cross-section. 
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9.7 Simulating and analyzing hadron-level events 

We are now in place to generate full events and analyze them with Sherpa. For simplicity 
we will not include a detector simulation in this exercise. However, Sherpa can be com- 
bined e.g. with PGS to simulate detector effects. For details on this procedure, please refer 
to the online manual, section 8.5 

http://sherpa.hepforge.0rg/doc/SHERPA-MC-l.4.O.html#PGS-interface. Sherpa also has 
a built-in Rivet-interface and it can output events in HepMC format. For more details, 
please refer to the online manual. 

Change back to the original directory and then go to ToyModel_HadronLevel. Have a 
look at the Run.dat input file. 

In the (run) section we have removed the switches that disabled parton showers and 
fragmentation. Instead there is a new switch, instructing Sherpa to perform a simple 
analysis. 

Open the Analysis.dat input file. It contains instructions for the built-in analysis 
module. 

• Finder 93 20 -4.5 4.5 0.4 1 

Construct kr-jets (kf-code 93) with D = 0.4, pT > 20 GeV and |r/| < 4.5. 

• Finder 11 15 -2.5 2.5 

Reconstruct electrons (kf-code 11) with pT > 15 GeV and \r]\ < 2.5. 

• DRMin 11 -11 0.2 

Require electrons to be separated from each other by AR > 0.2. 

• DRMin 11 93 0.4 

Require electrons to be separated from jets by AR > 0.4. 

Finally, we analyze the di-jet invariant mass distribution in the range < rrijj < 2000 on a 
linear scale with 100 bins: 

Mass 93 93 2000 100 Lin LeptonsJets 

Run this setup using 

<sherpa_pref ix>/bin/Sherpa 

As we have linked the Process/ and Results/ directory from the previous run, Sherpa 
will immediately start generating events. 

Once it has finished, plot the results of the analysis using 

. /plot_results . sh 

and have a look at the invariant mass distributions in plots. ps. Page one shows the 
changes when going from a pure parton- level simulation to hadron level, while page two 
has separate analyses of the various sub-processes. 
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9.8 Computing the partial widths 

Change to the directory Widths and open the input file Run.dat. 
The (processes) section contains setups for the following decays 

• U ^ ucf)! 

• U ^ U(l>2 

• E ^ e~^i 

• (j)2 e~e'^(j)i 

Comix is activated in the (me) section using COMIX_ALLOW_BSM 1. 
Let Sherpa compute the widths using 

<sherpa_pr ef ix>/bin/Sherpa 

Update the widths in the file Particle.dat. This file is linked from . ./FeynRules_Output 
and will be used in all toy model setups. 

9.9 Z+j backgrounds at LO and at NLO 

We will now proceed to simulate some important Standard-Model backgrounds. Sherpa 
will run substantially longer than in the previous steps. You may consider doing the fol- 
lowing part of the tutorial in parallel with another tutorial or after the workshop. 

Change back to the original directory and then go to Backgro\mds_ZJets. Have a look 
at the Run.dat input file. 

In the (processes) section you find the following: 

Process 93 93 -> 11 -11 93{2} 
Order_EW 2; CKKW sqr (30/E_CMS) ; 
End process; 

These lines instruct Sherpa to generate the process pp — t- e+e~ with up to two additional 
light partons. The tag CKKW indicates that the various sub-processes are to be merged using 
the truncated parton shower scheme. sqr(30/E_CMS) sets the value of Qcut to 30 GeV. 
Run this setup using 

<sherpa_pref ix>/bin/Sherpa 

Once Sherpa has finished, plot the results of the analysis using 

./plot_results . sh 

and have a look at the invariant mass distributions in plots. ps. 

If you like, check out the difference between a tree-level prediction of the background 
and the respective MENLOPS result. Run the MENLOPS setup using 

<sherpa_pref ix>/bin/Sherpa -f Run.NLO.dat 
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Once Sherpa has finished, plot the results of the analysis again, using 
. /plot_results . sh 

and have a look at the invariant mass distributions in plots . ps. Note that the NLO predic- 
tion is smoother because we have generated enhanced weighted events, cf. the Run . NLO . dat 
file. 

9.10 Top backgrounds 

Change back to the original directory and then go to Backgrounds_TTBar. Have a look at 
the Run. dat input file. 
Run the setup using 

<sherpa_pref ix>/bin/Sherpa 

Once Sherpa has finished, plot the results of the analysis using 

. /plot_results . sh 

and have a look at the invariant mass distributions in plots. ps. 

Now you can devise a strategy to reduce the Standard-Model backgrounds. 

A The history of the MC4BSM series of workshops 

The aim of the MC4BSM series of mini-workshops was to bring together theorists and 
experimentalists who are interested in developing and using Monte Carlo tools for BSM 
Physics in preparation for the analysis of data from the LHC. Since simulation of low energy 
supersymmetry (and the MSSM in particular) was already available from several excellent 
tools, the focus of the MC4BSM workshops (as the name suggests) was on simulation tools 
for alternative TeV-scale physics models. The main goals of the workshops were: 

• To survey what is available and provide user feedback on experiences with Monte 
Carlo tools for BSM. 

• To identify promising models (or processes) for which the relevant tools are still 
missing and start filling up those gaps. 

• To propose ways to streamline the process of going from from theory models to col- 
lider events (see Fig. 1). If this process were made more user-friendly, then a larger 
community of people would become involved in serious collider studies beyond the 
MSSM. 

The previous MC4BSM workshops were held at Fermilab (2006), Princeton (2007), 
CERN (2008), UC Davis (2009) and Copenhagen (2010). Many new tools (e.g. Bridge 
[41], Marmoset [42] Catfish [43], Pythia-UED [44], etc.) and model implementations 
(e.g. the left-right twin Higgs model [45, 46], the Littlest Higgs model with T-parity [47], 
Minimal UED [48], etc.) have been reported at these workshops. The latest, 6th installment 
in the MC4BSM series took place on the beautiful Cornell University campus in Ithaca, 
NY in March 2012. The weather was unseasonably warm and contributed to the lively 
atmosphere at the workshop. 
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